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We have measured the drift velocity (vp) of 
positive ions in liquid helium II in the tempera- 
ture range 1.06 to 0.927°K as a function of the ap- 
plied electric field. We have found that above a 
critical field strength E,. the drift velocity no 
longer increases linearly with the increasing 
field, but suffers a discontinuity. While the value 
of this critical field E,, depends on the tempera- 
ture, the corresponding critical velocity (v_) ap- 
pears to be a temperature -independent quantity 
in the range so far investigated. 

The apparatus used to measure the mobility is 
a time-of-flight method, which has been more 
fully described by one of us.’ Essentially one 
measures the drift velocity from the time taken 
by the ions to go from a grid to an electrode, 
which is connected to a vibrating reed electrom- 
eter by a suitable filter. Two apparatus have 


been used, with different electrode distances 
and different ionization sources, and they gave 
identical results. 

The reproducibility of our drift velocity data 
was within 2%, while the mobility up =(vp)/E 
derived in this way is affected by a further com- 
mon systematic error, due to the uncertainty in 
the absolute value of the square wave voltage 
which was estimated to rise up to 2%. In this 
range our absolute values of the zero-field mo- 
bility were reproducible on different days, and 
were also in agreement with those measured by 
Reif and Meyer’ and by ourselves® by different 
techniques. 

Some typical results of this investigation are 
reported in Fig. 1 in the customary way, as a 
plot of the mobility versus the applied field. A 
glance at this graph shows that below the critical 
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field of 37 v/cm the mobility is indeed field inde- 
pendent, but above it the mobility sharply de- 
creases to a lower constant value. The sharp 
change in the mobility is well outside the experi- 
mental error, and also the discontinuity is clear 
from the shape of the curve around E-. 

In Fig. 2 we give the results of some runs made 
at different temperatures in a plot of mobility 
versus drift velocity, particularly convenient to 
show that the discontinuity always falls around 
the same value of the drift velocity, while the 
mobility changes by a factor of 3. For the sake 
of clarity the experimental errors in the mobility 
and drift velocity are not shown in this plot. In 
Table I we give the critical velocities so far de- 
tected under different experimental conditions. 
Each value of (v¢) is evaluated as the mean value 
of the two measured drift velocities where the 
break was found; the uncertainty in (v c) indi- 
cated in the table is therefore one half of the 


Table I. Observed critical velocities. 











Run T (°K) (v) (m/sec) 
3B 1.036 4.79 + 0.10 
3A 1.000 4.83 +0.10 
5A 0.966 4.56 +0.15 
4A 0.928 4.77+0.10 
8A 0.927 4.56 + 0.10 





value of this range. These values of (v_) cluster 
in the range (v,) =4.70+ 0.15 m/sec, indicating 
indeed the existence of a threshold independent 
of the temperature. 

In one run we were able to achieve still higher 
drift velocities, and a second discontinuity was 
observed as shown in Fig. 3. This second dis- 
continuity appeared at a value (v-) = 9.22+ 0.20, 
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FIG, 3. A similar plot 
to Fig. 2 for one temper- 
ature and for an extended 
range of drift velocities. 
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which is larger than the first one by a factor of 
1.96+0.08, as if a repetition of this phenomenon 
can exist at a velocity twice as large. 

The change in the mobility values at the discon- 
tinuity is never larger than 10% and therefore 
this phenomenon may easily escape observation 
in a plot of the mobility versus field when few 
experimental points are taken, and may then be 
confused with a slight and continuous field de- 
pendence. We believe this to be the reason why 
Reif and Meyer’ did not notice this phenomenon 
at 0.98°K, where they were interested only in the 
zero-field mobility. It is also clear that at tem- 
peratures higher than the ones reported here it 
becomes progressively harder or downright im- 
possible to detect this discontinuity. Some pre- 
liminary runs performed in this laboratory, at 
temperatures lower than 0.90°K, indicate that 
this discontinuity in the mobility is more diffi- 
cult to detect there because of the superimposed 
continuous decrease. Reif and Meyer” have con- 
centrated their work on the field dependence of 
the mobility essentially below this temperature, 
and have thoroughly investigated this continuous 
decrease in terms of mean-free-path phenomena. 
Therefore we believe our experimental work not 
to be in contrast with their data, but instead to 
provide a fine structure of the mobility versus 
field curve which is observable only in a narrow 
temperature range. 

An interpretation of the meaning of the critical 
drift velocities detected in these experiments is 
not easily found in the current picture‘ of the ionic 
motion in liquid helium Il. The sharpness of the 
transition and the stepwise shape of the curve 
rules out any approach in terms of hot ions, which 
is customary to explain the decrease of the mobil- 
ity with the increasing field strength when the 
drift velocity is comparable to the thermal veloc- 
ity of the ions. On the other hand, we know that 
these velocities are much too low to be possibly 
explained in terms of inelastic processes involv- 
ing the creation of one roton or phonon. 

As a further possibility one may suggest the 
same argument now used to explain the critical 
flow phenomena, which involves the creation of 


quantized vortex lines. The creation of one quan- 
tized vortex ring of atomic size behind every ion 
in motion would easily provide a critical velocity 
of the magnitude observed here, and would also 
give an account of the observed multiplicity in 
terms of different quantum numbers of circula- 
tion. However, if the Feynman® expression for 
the vortex ring energy and momentum are used, 
the effective mass of the ion turns out to be as 
large as 10* helium atom masses. Details of the 
calculation are not given here, but this figure 
follows simply from the energy and momentum 
conservation conditions.” Such a value of the ef- 
fective mass is much larger than the estimated*’’ 
static mass, and seems rather inconsistent. 

One final remark should be made concerning 
the unexpectedly sharp onset of the drop in mo- 
bility. The only two possible explanations which 
suggest themselves are (1) a very sharply peaked 
velocity distribution due to a large effective mass 
of the ion, or (2) a kind of cooperative interaction 
with a large number of nearby ions through the 
long-range Coulomb forces. ; 

The lack of a proper theoretical treatment of 
this problem does not allow us at the present time 
to come to a more definite conclusion. 

Experiments are planned in this laboratory to 
extend these measurements over a wider range 
of field strengths and temperatures by an im- 
proved technique. 





1§, Cunsolo, Nuovo cimento (to be published). 

°F. Reif and L. Meyer, Phys. Rev. 119, 1164 (1960). 

5G. Careri, F. Scaramuzzi, and J. O. Thomson, Nuo- 
vo cimento 13, 186 (1959). 

‘For a discussion on the behavior of ions in liquid he- 
lium, see also G, Careri, Progress in Low-Tempera- 
ture Physics (North-Holland Publishing Company, Am- 
sterdam, 1961), Vol. 3, p. 58. 

5k. P. Feynman, Progress in Low-Temperature 
Physics (North-Holland Publishing Company, Amster- 
dam, 1955), Vol. 1, p. 17. 

*k, R. Atkins, Phys. Rev. 116, 1339 (1959). 

"J. de Boer and A. ‘t Hooft, Proceedings of the Sev- 
enth International Conference on Low-Temperature Phys- 
ics (University of Toronto Press, Toronto, 1960), p. 
510. 
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CALCULATED SHAPE OF THE FERMI SURFACE OF COPPER 


B. Segall 
General Electric Research Laboratory, Schenectady, New York 
(Received July 12, 1961) 


The electronic properties of copper have recent- 
ly become the subject of intensive study by a wide 
variety of experimental methods.’ One important 
result of these investigations has been the general 
verification of the Pippard model* for the Fermi 
surface, which is nearly spherical over the main 
part of the surface (called the “belly”) and which 
is pulled out in the (111) directions to form eight 
“necks” which contact the Brillouin zone surface. 
It is the purpose of this Letter to note that our 
recently completed band-structure calculations*® 
on Cu not only have led to results consistent with 
the Pippard model but have also predicted interest- 
ing distortions of the “belly” which have just been 
observed experimentally.* Also we will present 
an explanation for the distortions and will discuss 
the importance of this effect for the other noble 
metals. 

The calculations were carried out by the Green’s 
function method’ using two different crystal poten- 
tials. The first one, which is perhaps most ap- 
propriate for the d bands, was used initially by 
Chodorow.® For the second potential, we com- 
puted the Coulomb and exchange contributions 
separately for the different angular momenta from 
renormalized Cu Hartree-Fock functions. Only 
the results for the latter potential, which we be- 
lieve to be somewhat more realistic than the for- 
mer, will be discussed here although the other 
results are in semiquantitative agreement with 
them and lead to similar conclusions about the 
Fermi surface.’ 

Our results® show that the Fermi surface as- 
sociated with the calculated energy bands contacts 
the hexagonal zone surface and the calculated ra- 
tio of neck to belly radii is 0.22+ 0.02 in good 
agreement with the experimental ratio’»* of about 
0.20. Also, for these bands, the onset of inter- 
band optical transitions would occur for transi- 
tions between the d bands around the upper L, 
state and the Fermi surface at about 2.6 ev in 
fair accord with the observed 2.2-ev value. 

To study the Fermi surface in more detail, we 
have determined the intersections of the constant 
energy surfaces with a (110) plane (Fig. 1). The 
estimated Fermi energy is -0.19 ry. Aside from 
the sizable region of contact with the zone sur- 
face (for E=E F)s the principal feature of interest 
in these curves is that the belly is pulled out in 
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the (100) directions and is contracted in the (110) 
directions. The occurrence of these appreciable 
distortions is at variance with the generaily held 
belief that the main section of the surface is very 
nearly spherical. 

Bearing out these results are the recently re- 
ported studies of the shape of the Cu Fermi sur- 
face by Morse et al.* using the magnetoacoustic 
effect which clearly show these deviations from 
sphericity. It is clear from their diagram that 
their observed and our calculated distortions are 
comparable. We find k -(100)/k (110) = 1.09. 

The explanation for the occurrence of these dis- 
tortions involves the “repulsion” of conduction 
band states by the relatively high-lying d-band 
states of the same symmetry. In Fig. 2 are shown 
some of the bands of interest for this discussion. 
Now, the calculations® show that the lowest con- 
duction band state at the end of the [100] axis has 
the p-like symmetry X,’. As d-band states of this 
symmetry cannot occur, this state is not “re- 








E=—,187 ry 








FIG, 1. The intersections of a (110) plane with cal- 
culated constant energy surfaces for Cu and with a free 
electron Fermi sphere (dashed curve). The estimated 
Fermi energy is -0.19 ry. 
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FIG, 2. The Cu energy bands in the conduction and 
d-band regions for the A;, 4;, and 2, states in the 
[111], [100], and [110] directions, respectively. For 
clarity, most of the calculated d-band curves were 
not drawn in. 


pelled.” Conduction band states “near” X,’ will, 
of course, have very small admixtures of d func- 
tions and consequently will only be affected weak- 
ly by the lower states. In the [110] direction, on 
the other hand, the conduction band (=, symmetry) 
“interacts” fully with the d band all the way out to 
the zone surface where it has K, symmetry. The 
argument also applies to the other noble metals 


as our calculations for them® also indicate that 
the lowest conduction band states at X and L are 
p-like. It is interesting to note that Morse et al. 
have also observed these distortions in Au. Inas- 
much as the d bands in Ag are expected to be low- 
er than in Cu and Au, we would predict smaller 
distortions in that metal. 

It is thus seen that, contrary to opinions often 
expressed in the past about the possibility of 
calculating meaningful energy bands for such a 
metal, the calculated bands lead to a Fermi sur- 
face in good agreement with the one determined 
by experiment. Of course, as the many-elec- 
tron effects have not been treated adequately and 
the electron-phonon coupling has been neglected, 
we cannot expect these bands to describe all the 
electronic properties of this metal without the in- 
troduction of corrections for these correlation 
effects. 

The details of the calculations will be published 
shortly. 

The author wishes to thank Miss E. L. Kreiger 
for her invaluable help with the numerical compu- 
tations. 





'This work is reviewed in The Fermi Surface, edited 
by W. A. Harrison and M. B. Webb (John Wiley & Sons, 
Inc., New York, 1960), in the articles by D. Shoenberg, 
p. 74; J. R. Klauder and J. E. Kunzler, p. 125; A. F. 
Kip, p. 146; G. E. Smith, p. 172; and R. W. Morse, 
p. 214. 

2A. B. Pippard, Phil. Trans. Roy. Soc. London 
A250, 325 (1957). 

The essential results and conclusions discussed in 
this paper were given in B. Segall and E. L. Kreiger, 
Bull. Am. Phys. Soc. 6, 10 (1961), and in B. Segall, 
Bull, Am. Phys. Soc, 6, 231 (1961). The latter talk 
has been written up as General Electric Research Lab- 
oratory Report No. RL-2785G which is available upon 
request. Also, see B. Segall, reference 1, p. 314. 

4R. W. Morse, A. Myers, and C. T. Walker, J. 
Acoust. Soc. Am. 33, 699 (1961). 

‘w. Kohn and N. Rostoker, Phys. Rev. 94, 1111 
(1954), and F, S. Ham and B. Segall (to be published). 

81, Chodorow, thesis, Massachusetts Institute of 
Technology, 1939 (unpublished). 

Tore recently G. A. Burdick (private communica- 
tion) has aiso calculated E(k) for Chodorow’s potential 
and his results agree with ours. 

8B, Segall, Bull. Am. Phys. Soc. 6, 145 (1961); 6, 
231 (1961). 
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TOPOLOGY OF THE FERMI SURFACE OF COPPER 


Glenn A. Burdick 
Massachusetts Institue of Technology, Cambridge, Massachusetts 
(Received July 7, 1961) 


Segall’ has reported calculations of the Fermi 
surface of copper by the Green’s function method 
of Kohn and Rostoker, using both the potential of 
Chodorow, ? and other modified potentials. In this 
laboratory® we have been engaged since 1954 in 
programming the augmented plane wave method 
suggested by Slater* in 1937 for the digital com- 
puter, and in making calculations of energy bands 
by this method. The author has used these pro- 
grams to compute the energy versus k& for copper, 
for the equivalent of 2048 points in the first Bril- 
louin zone and for the highest 6 bands which are 
fully or partially filled. The potential derived by 
Chodorow for d electrons was used in the compu- 
tations, which were done on an IBM-709 computer 
at the M.I.T. Computation Center. 

The purpose of this report is to point out that we 
are in agreement with the results of Segall, name- 
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ly, that there is very close agreement between the 
Fermi surface for copper as proposed by Pippard® 
and the one as theoretically computed by the author. 
The major features of this surface are the “necks” 
which project out from the “belly” region in the 
(111) directions and intersect the Brillouin zone 
edge at the center of the hexagonal faces. The 
author also finds a slight bulge (in the form of an 
“egghead”) in the (100) directions. 

Figure 1 shows the Fermi surface cross section 
as formed by a plane passing through the origin 
and bisecting two opposite sides of a hexagonal 
face. Several neighboring zones have been drawn 
to facilitate the definition of characteristic dimen- 
sions of the surface. The experimental values for 
the dimensions as defined in Fig. 1 were obtained 
from ultrasonic attenuation experiments® and are 
listed along with the author’s theoretical values in 


(110> > 


FIG, 1. The (110) cross sec- 
tion for the Fermi surface of 


copper. 
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Table I. Experimental values® for the dimensions 
as defined in Fig. 1. For comparison, the radius of 
the sphere having one-half the volume of the Brillouin 
zone is k, =1.44. 











Morse Burdick 
Ry00 1.41 1,.5120.03 
R410 1.38 1,.39+0.03 
kq 1,22 1,21+20.03 
kn 1.10 1.1120,.02 
ke 0.32 0.44 +0.06 
ky 0.24 0.24+0.03; 0.22+0.06 








*all values are to be multiplied by the factor 107". 


Table I. The possible computational errors are as 
indicated but the author has no information concern- 


‘ing the possible experimental errors. The dimen- 


sions are given in cgs momentum units. The two 
theoretical values listed for k, correspond to the 
“neck” radius as computed along a line passing 
through the center of the hexagonal face and one 
corner and along a line passing through the center 
and the midpoint of one edge, respectively. As 
can be seen from an examination of Table I, the 
only difference in topology is that the author’s 
“belly” region in the plane of the figure is slight- 
ly more egg-shaped than the experimental model 
indicates. However, even this difference is a 


slight one. 

As those who heard the two talks of Segall before 
the American Physical Society (referred to above) 
will realize, these results agree closely with his, 
thereby giving important verification of the fact 
that energy bands as computed by different meth- 
ods can agree with each other and with experiment 
to a greater extent than has been generally rea- 
lized. 

The details of the computations, the density-of- 
states curve, and a discussion of the energy band 
structure will be published at a later date. 





'B. Segall and E. L. Kreiger, Bull. Am. Phys. Soc. 
6, 10 (1961); B. Segall, Bull. Am. Phys. Soc. 6, 231 
(1961). We are greatly indebted to Dr. Segall for pri- 
vate communication of many of his results. 

2M. I. Chodorow, Ph.D. thesis, Massachusetts In- 
stitute of Technology, 1939 (unpublished). 

SQuarterly Progress Reports, Solid-State and Mo- 
lecular Theory Group, April 15, 1954 to date (unpub- 
lished), contributions from D. J. Howarth, M. M. 
Saffren, J. H. Wood, and others. 

‘J, C. Slater, Phys. Rev. 51, 846 (1937). 

°a, B. Pippard, Phil. Trans. Roy. Soc. London 
A250, 325 (1957). 

Robert W. Morse, in The Fermi Surface of Metals, 
edited by W. A. Harrison and M. B. Webb (John Wiley 
& Sons, Inc., New York, 1960), p. 214 ff. 








CALCULATION OF THE CRYSTAL FIELD SPLITTING 


R. G. Shulman and S. Sugano 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received July 21, 1961) 


By starting with wave functions determined by 
nuclear magnetic resonance measurements? in 
KNiF,, we have been able to calculate a value 
for the cubic crystal field splitting which is in 
good agreement with the observed value. The 
cubic crystal field splitting, 10Dqg, is the energy 
separation between the orbital triplet to, and the 
€g doublet which arise from the d levels split by 
a cubic field. Measured values of 10Dq have 
been used with considerable success to explain 
Optical and magnetic properties. However, it 
has not yet been possible to calculate a value 
which agrees with experiment. We follow the 
method of Tanabe and Sugano,” the main differ - 
ence being the use of molecular orbitals with 
Hartree-Fock atomic functions as our basis set. 






We define 
s 5,33 a 6,23 
10Dq E(t, € T,) E(t, e A,) 
= h a 1 
(wv! Iw) ( lnlv), (1) 
in which the wave functions are assumed to be 
=N “V2 o 
°, e [e, r Yeh (2) 
=N “v2 - q 
In these expressions $, and 4, (subscript ¢ is ab- 
breviation for ¢,) may be the Ni** 3d,2_,2 and 
3dyy functions, respectively, with Watson’s Har- 


tree-Fock radial functions.* Furthermore, ¥, 
=UgeXo+HseXs and Uge"+Hse*=1 while xq; Xs; 
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Table I. Values of the matrix elements listed in the 
first column are shown in the second in atomic units. 
In the third column the contributions to 10Dq from the 
differences between the eg and tog qnenyies are listed 
incm=!, The calculated value of 5270 cm™! agrees 


fairly well with the experimental® value of 7250 cm™, 





Contribution to 





Matrix element Value (a.u.) 10Dq (em™') 
(® |al@ ) +0.1201 
e e -3580 
(@ |nl®,) +0.1364 
(@ |hlv ) -0.1510 
é é +16 320 
(@ |al¥,) -0.0559 
(Y |hlv ) -0.3300 
Oe. -7470 
(¥ lnl¥) -0.0847 
Calculated 10Dq +5270, 
Observed 10Dq +7250 





% | Knox, R. G. Shulman, and S. Sugano, Bull. Am. 
Phys. Soc. 5, 415 (1960). 


and x¢= yz are the normalized appropriate linear 
combinations of fluoride ion 2p,, 2s, and 2p, or- 
bitals, respectively. The fluoride ion 2s and 2p 
atomic orbitals were Froese’s Hartree-Fock func- 
tions.* The fluorine 1s functions were shrunk into 
the nuclei. The augmentation parameters for ¥, 
Of AgU ge =Ag and Agi se =Ag were determined from 
the nuclear magnetic resonance measurements’ 
of the F’* nucleus in a single crystal of KNiF,. 
The values® were (Agi ge)” =3 x 0.0369 and (Ag se)? 
=3 x 0.00451. 

The first order Hartree-Fock Hamiltonian em- 
ployed in (1) is 


h=-44+V (4) 


V 
ML’ 
in which A is the kinetic energy operator and Vjy 
is the Coulomb and exchange potential of all the 
Ni** electrons and nucleus. The contribution of 
the ligands is 


0 K E 
v,7V, +V, +V, ; (5) 
in which 
6 
*s “2 irk, 1'7~ eT 


is the point-charge Pn of six ligands cen- 
tered at Rj; 


6 
K 8 -1 2 
Vv. = “re | Yr >» 14, (2)! 
L i=l ly RI 2 12 k=2s, 2p ik 

is Kleiner’s correction® of the point-charge poten- 
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Table II. Separate contributions to the nickel diagonal 
terms of the components of V; . 








0 K E 
ve ve ", 
(@ |v_|@) 1.58687 -0.03017 -0.02396 
e Le 
(@ |v, 1%) 1.58055 -0.02069 -0.01082 
Difference 0.00632 -0.00948 -0.01314 








tial which arises from the imperfect screening of 
the fluorine nuclei by their electrons whose atom- 
ic orbitals are $;,; finally 


= fr : | qr. io" 2%, °(2)%,, (1)P 


k=2s, 2p 12 


is the exchange interaction introduced in the pre- 
vious paper,” where P,, is the permutation opera- 
tor of electrons 1 and 2. All the two-center inte- 
grals were calculated on an IBM-7090 with Swit- 

endick and Corbat6’s program.” The results are 

presented in Table I. 

The nickel diagonal term makes a negative con- 
tribution to 10Dq. This term represents the ap- 
proximation in which the effects of the ligands 
are treated as perturbations upon the nickel ion 
energy levels. The separate contributions to 
this term of v,°, VX, and v p= are shown in 
Table II. Although the individual exchange terms 
are small, they make a larger contribution than 
the total Coulomb interaction because of the can- 
cellation® of the contributions from V 1° and V; 
which allows the difference between exchange 
terms for e and ft, to dominate. In Table II, the 
difference between N; and N¢ is not taken into ac- 
count because its effect is negligible. 

In order to calculate the contributions of the 
transfer integral and fluorine diagonal terms to 
10Dq, it is necessary to know the value of A; 
since the nuclear magnetic resonance measure- 
ments only determine A, andA,. Az may be de- 
termined nonempirically by minimizing the cor- 
responding bonding orbital energy. Such a non- 
empirical method can be also applied to the de- 
termination of A, andA,, and a comparison be- 
tween their calculated and observed values would 
be very interesting. 

In the present calculation, however, we em- 
ployed the following method in order to avoid the 
tedious calculation of important three-center in- 
tegrals of the type ($3q1Vjz1;z) appearing in 
the transfer integrals; a three-center integral 
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for the e orbital was determined so that the ob- 
served A, minimizes the e bonding orbital ener- 
gy. Then, using the calculated overlap integrals 
to relate the three-center integral for the ¢, or- 
bital to that for the e orbital, A; was determined 
by minimizing the energy, with the result that 

dy =0.286. The large value of Az is consistent 
with our previous findings.’ Furthermore, the 
values of the three-center integrals thus deter- 
mined were found to be reasonable. The three- 
center integrals are being evaluated independently 
so that we will be able to compare calculated and 
observed values of Ag andA,s. However, the rea- 
sonable values derived for the three-center in- 
tegrals in the above-mentioned procedure con- 
firm the observed values of A, and A,. It is in- 
teresting to observe that the V Le term makes an 
important contribution to the 10Dgq in the trans- 
fer integral also. 

In all energy determinations we have assumed 
a potential well model in which the NiF,~* com- 
plex is located in a constant potential arising 
from the rest of the crystal. This model can be 
proven to be a good approximation. It also ex- 
plains the observation that the cubic field param- 
eter is insensitive to distant ions beyond the lig- 
ands. All of these assumptions, calculations, 
and additional refinements will be discussed 
more fully in future publications. 

In conclusion we can say: 

(1) The traditional concept of the cubic crystal 
field should be changed: The nickel diagonal 
term gives a negative contribution to the crystal 
field splitting. Its contribution comes mainly 
from the exchange interaction between d electrons 
and ligand electrons. The transfer integral term 
is the only term which makes a positive contribu- 
tion to 10Dq. Roughly speaking, this contribution 
is reduced by a factor of ~50% by the negative 


contribution of the fluorine diagonal term. Al- 
though this relation is found in our specific prob- 
lem numerically, it is generally applicable. This 
will be discussed in the future. 

(2) The molecular orbital model is successful 
in explaining both the nuclear magnetic resonance 
experiments and the cubic crystal field param- 
eter in a consistent way. The value of 10Dq cal- 
culated from this model is 5300 cm™, in reason- 
able agreement with the observed value of 7250 
em™, 

It is a pleasure to thank A. C. Switendick, 
Mrs. W. Mammel, and R. Kornegay for help 
in the machine calculation and L. C. Allen for 
an important suggestion. 
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FERROMAGNETIC INTERACTION IN EuO 


B. T. Matthias and R. M. Bozorth* 
Bell Telephone Laboratories, Murray Hill, New Jersey 


and 


J. H. Van Vieck 
Harvard University, Cambridge, Massachusetts 
(Received August 3, 1961) 


We have found that a europium oxide of the for- 
mula EuO becomes truly ferromagnetic at 77°K 
with a saturation moment of close to 7 Bohr mag- 
netons. This is thus the first rare earth oxide 
to be found to become ferromagnetic, and with 
the exception of CrO, the only oxide to our knowl- 
edge that has true ferromagnetic coupling. 

We were led to this discovery by somewhat con- 
tradictory results on ferromagnetic intermetallic 
rare earth compounds.’ There it was found that 
Eulr, becomes ferromagnetic near 70°K and has 
an extremely small saturation moment. It was 
clear from measurements of lattice constants of 
this cubic compound that europium had a 4f°* con- 
figuration corresponding to trivalent europium and 
so would not be expected to have a ferromagnetic 
moment.* A compound having a Curie point of 70°K 
and very small moment seemed to be a severe 
contradiction. One possible explanation was the 
presence of an impurity strongly ferromagnetic 
below 70°K. Eulr, had been formed in a quartz 
tube by having molten europium react with fine- 
ly divided iridium powder, and the presence of ox- 
ides was thus a likely possibility. The magnetic 
behavior of Eu,O, is known and normal. EuO was 
the other alternative. 

The latter compound was formed by reacting 
stoichiometric amounts of Eu and Eu,O, at 1220°C 
in a quartz tube. The product was a brown crys- 
talline powder. X-ray examination by Mrs. V. B. 
Compton showed it to have a NaCl structure with a 
lattice constant of 5.141 A, corresponding to a den- 
sity of 8.21. A few unidentified faint lines were 
present on the x-ray film. This lattice constant 
agrees with that previously reported for material 
obtained by the reduction of Eu,O, with lanthanum 
metal.°® 

The existence of strong ferromagnetism in EuO, 
and a Curie point (77°) very near to that observed 
for Eulr, (70°) shows immediately that the presence 
of a few percent of EuO might easily give rise to 
the anomalous results obtained with Eulr,. 

As shown in Fig. 1, above the Curie point the sus- 
ceptibility of EuO is represented within the experi- 
mental error by a formula of the type x=C/(T-4), 
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where the Curie constant C has the value 7.9N§?/3k 
to be expected for an °S configuration. This fact 
shows unequivocally that the europium ion is diva- 
lent and isoelectronic with Gd*** as one would 
expect from the degree of oxidation. The line- 
arity of 1/x in T is, if anything, more accurate 
than one would expect, for the Weiss-Curie law 
is always but an approximation theoretically. 
Figure 2 shows the data for the ferromagnetic 
region. The saturation at 0°K, determined by 
extrapolation linear in 1/H, is 6.9 Bohr magnetons, 
and within the limits of error is equal to the theo- 
retical Jg=7. The broken line is that calculated 
for spontaneous magnetization using the Brillouin 
function for J=7/2 with the parameter a= (7/3) 
x (M/M,)/(T/@). At low temperatures the triangu- 
lar points are obtained by extrapolation (using 1/H) 
of the smoothed data for H =8 and 12 koe at se- 
lected temperatures. Just below the Curie point 
(obtained by the Weiss-Forrer method, W.F. in 
the figure) the extrapolation is to H =0, on the as- 
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FIG, 1. Reciprocal molecular susceptibility of Eu0: 
experimental points, and theoretical line drawn through 
ferromagnetic Curie point (77°K) with slope for Jg=7, 
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FIG, 2. Magnetization of EuO vs temperature in 
ferromagnetic region. Broken line, Brillouin curve 
for J= i, 6=77°. Triangles, extrapolation of experi- 
mental data to H== (high M) or H=0 (lowM). W.F., 
Weiss-Forrer extrapolation (¢,,) 7 = 9 Vs T, to deter- 
mine the Curie point. 


sumption that here the crystal anisotropy is zero 
and the applied field increases the magnetization 
above the spontaneous value. 

At low temperatures the magnetic moment per 
unit volume, based on 78 and a density of 8.25, 
is high: 

4nM = 24100. 
This may be compared with the values of 47M for 





iron (21700), gadolinium (about 25000), and for 
metallic dysprosium* (about 38 000). 

Measurement of the pressed powder on a meg- 
ohm meter indicated a resistivity of 10° to 10° 
ohms. 

The measurements are of interest as furnish- 
ing one of the clearest examples of exchange in- 
teraction in a material in which the atomic orbit- 
als are substantially the same as for the free 
ions, since this is always the case for 4f elec- 
trons in the middle of the rare earth period in 
nonconducting compounds. 

We will not attempt to enter into the question 
of the precise mechanism responsible for the 
ferromagnetic exchange coupling. It is certainly 
not the Yosida mechanism resulting from the 
polarization of conduction electrons, since EuO 
is an insulator. Conceivably, direct overlap may 
be the cause, but in our opinion indirect exchange 
via excited states is more likely. 

We are grateful to Mrs. V. B. Compton for per- 
mission to quote the results of her x-ray meas- 
urements on Eulr, and EuO. 
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QUANTIZATION OF FLUX IN A SUPERCONDUCTING CYLINDER 


John Bardeen 
University of Illinois, Urbana, Illinois 


(Received July 31, 


Recent experiments of Deaver and Fairbank’ 
and of Doll and Nabauer? have shown that the flux 
in a superconducting cylindrical tube is quantized. 
That the flux should be quantized in units of hc/e 
was originally suggested by London*; the experi- 
ments indicate units of about half this amount. 
Onsager* predicted such a result on the basis that 
an effective charge e*= -2e, representing that of 
a pair, should be used. We shall show that in 
tubes of very small diameter and with wall thick- 
ness of the order of the penetration depth the unit 
may depend on dimensions and be somewhat small- 
er than hc/2e. 

The theory has been discussed by Byers and 
Yang,° who showed that it is essential to take into 
account the statistical distribution of qudsi-parti- 
cle excitations, and that one could expect a flux 
quantum of hc/2e on the basis of the microscopic 
theory. They considered, as we shall do here, 
states with current flow which exist with nonin- 
tegral values of flux. 

We shall use the Gor’kov version® of the Ginz- 
burg-Landau (G-L) theory’ to discuss the effect. 
By use of the method of “thermal” Green’s func- 
tions, Gor’kov showed that the G- L equations, 
originally given on a phenomenological basis, may 
be derived from the microscopic theory provided 
that the temperature is sufficiently close to T, so 
that the local London equations are valid for the 
electrodynamics. The only difference from the 
original version is that an effective charge e* 
=-2e, evidently that of a pair, appears in place 
of -e. Since in the experiments, the flux is fro- 
zen in as the specimen is cooled below T,, the 
theory is valid in the significant temperature 
range. 

Gor’kov showed that the effective wave function, 
v(r), of the G-L theory may be taken to be pro- 
portional to the local value of the energy gap. We 
shall take a different normalization for © than that 
used by Gor’kov, one that is closer to the original 
G-L theory. We shall define it so that ¥,, the 
equilibrium value of ¥, is such that |¥9|?=p,/p, 
where p, is the density of superconducting elec- 
trons of the two-fluid model and p=Nm, the total 
density. The Schrédinger-like equation for ¥ may 
then be written 


2 
“sae arte 4 arie4@)) ¥ - a¥+BIvl*¥=0, (1) 
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where m*=2m and e*=-2e, representing the 
charge and mass of a bound pair. This equation 
then corresponds to that of the center-of-mass 
motion of the pair. The free-energy difference 
between superconducting and normal states is giv- 
en by integrating over space and multiplying by 
N/2, the number of pairs at T=0°K: 


N h?__+/(8 2nie* ~\" 
F, -F,,-3f |-sae (e- he AG@) 





ali? $ol¥ Ian, (2) 
The density of mass flow is 
A ee er yl? 
J= i -#[v oY c|-< wra@l. (3) 


In the absence of fields and currents, the free- 
energy minimum, for |¥9|?=a/8=p,/p, is 

Fy, -Fy= -H,?/81 = -Na?/(48), where H-(T) is the 
critical field for bulk material. 

Gor’kov derived the values of a and f from the 
microscopic theory.® They can be expressed in 
terms of H, and the penetration depth, »(7). Since 
ps/p=Az?/a?, where \,7=mc?/(4ne7N) is the 
London penetration depth, we have 


a =H *r*/(20Nd,*); B=Ha*/(20NX,*). (4) 


These apply whether or not elastic scattering by 
impurities or other imperfections is important 
provided that the appropriate value for the pene- 
tration depth is used. Limiting expressions®: ° 
valid near T, are 


- os ‘ 5 
d/*/re= 2(1 t), l> E03 (5) 
d,7/r?= 2.73(1- t)(l/Eo), l«é); (6) 


where ¢t=7/T,, 1 is the mean free path, and £, is 
the coherence distance. The limit /« é, applies 
to thin films. 

In the absence of a magnetic field, a wave func- 
tion v=, eig-F corresponds to a momentum 
m*vs =hq/2 7 where vg is the velocity of the cen- 
ter of mass. The net current is p,v, and the in- 
crease in free energy 3p,v,”, in accordance with 
the two-fluid model.*° A superconductor is char- 
acterized by a value of p, different from zero. 
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To discuss the quantization of flux, we shall 
suppose that the specimen is in the form of a 
thin film of thickness d<2) on the surface of a 
cylindrical tube of radius y and circumference 
L=2ar. The current density will then be sub- 
stantially uniform across the film. One may sim- 
plify the mathematics by taking x to be a linear 
coordinate measuring distance around the cir- 
cumference and require periodic boundary condi- 
tions, ¥(x+L)=(x). If @ is the flux through the 
tube, the vector potential A, is a constant equal 
to 6/L. 

An appropriate periodic solution of the G-L 
equation is*® 


W = exp[27inx/LW,, (7) 


where n is an integer. For an arbitrary flux @ 
one may choose v so that the residual wave vec- 
tor 


q= (27/L)(n -e*&/he) (8) 


is in the interval -1/L<q<1/L. The velocity 
v, =hq/(2mm*) and the free-energy increase will 
then be as small as possible. 

The electric current density, I=e*p,v,/m*, is 
quite large even for small values of g. A value 
q=1/L, with L~1 cm, gives v,~1 cm/sec and a 
current density J~10° p,/p amp/cm? in a typical 
metal. Generally a much smaller current density 
is required to produce a unit of flux, so that it is 
favorable to minimize the free energy by reduc- 
ing g and thus vs to nearly zero. There will then 
be very nearly an integral number of flux units in 
the cylinder. Values of g corresponding to odd 
multiples of 7/L would also give zero current 
(with a different pairing condition), but with a 
free-energy maximum rather than a minimum.® 

Appreciable values of g which would give non- 
integral units of flux can be expected only if the 
thickness of the film and the radius of the cylin- 
der are very small. When the external field is 
removed, the flux is produced entirely by the cur- 
rent flowing in the film, so that 


$= 41?y7I1d/c. (9) 





Setting 1=e*p,hq/(81m?), we may solve for g and 
find 


q= 2,7e*@/(r*dhc), (10) 


where )?=m?*c?/(4me"p,). Substituting in (8) and 
solving for @, we find 


& = (nhc /2e)(1 + 2n?/rd)™. (11) 


Thus we find that for tubes of very small radius, 
the flux quantum may be appreciably less than 
hc/2e. It is possible that the value of 0.4hc/e ob- 
served by Doll and Nabauer may be accounted for 
in this way. The lead tube which they used had a 
diameter of only 0.01 mm and \ may have been 
abnormally large because of scattering in the film. 
It should be noted that ) and thus the unit of flux 
may change with temperature. 

This work was done at the Summer Institute for 
Theoretical Physics, University of Wisconsin, 
sponsored by the National Science Foundation. The 
author is indebted to E. Wigner, R. G. Sachs, and 
other members of the Institute for informative dis- 
cussions. 
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QUANTIZATION OF THE FLUXOID IN SUPERCONDUCTIVITY 


J. B. Keller and B. Zumino 
New York University, New York, New York 
(Received August 16, 1961) 


It was predicted by London’ that the magnetic 
fluxoid in a superconductor should be quantized 
in units of 2mfc/e. On the basis of the electron 
pair idea of Schafroth? and Cooper,® Onsager* 
suggested that the quantum unit should be half 
that of London, since the charge of a pair is twice 
the charge of the electron. Onsager’s prediction 
was confirmed in recent experiments of Deaver 
and Fairbank® and of Doll and Nabauer® on super- 
conducting cylindrical tubes. Theoretical discus- 
sions of the phenomenon have been given by Byers 
and Yang,’ Blatt,® and Bardeen.® Bardeen showed 
that the quantum of flux is smaller than the natu- 
ral unit tfic/e for a tube thinner than the penetra- 
tion depth. 

It is our purpose to show that in a superconduc- 
tor the fluxoid is always quantized in multiples of 
the natural unit even though the flux may be quan- 
tized in smaller units. Our demonstration is based 
upon the assumption that there exists an effective 
wave function © for electron pairs. We also as- 
sume that ||? is the number density of pairs and 
that the supercurrent density J is given by the us- 
ual quantum mechanical formula 


+ -ieh ,* o, 
a 7 (YW VY - VN poe A. (1) 


Here A is the vector potential and we have used 
the fact that the charge and mass of a pair are 
twice those of an electron. These assumptions 
have been verified by Gor’kov"® in his derivation 
of the Ginzburg-Landau"! semiphenomenological 
theory from the microscopic theory of Bardeen, 
Cooper, and Schrieffer.'” 

Following London, we define the fluxoid oo 
through a closed curve C by the equation 


0= [,(A+cas)-as, (2) 
where 
A=m/(2e?|¥|?). (3) 
From Eq. (1) it follows that 
A+cAJ= -(ihc /4e)¥V in(¥/¥"). (4) 
Then Eq. (2) becomes 
4, = -(iic/4e)[In(w/¥")], (5) 


In Eq. (5), square brackets denote the change in 
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the enclosed quantity upon traversing C once, 
From Eqs. (4) and (5) it follows that @> has the 
same value for any two curves which can be de- 
formed into each other without leaving the super- 
conductor. Therefore $c vanishes for a curve 
which can be shrunk to a point. For a curve which 
cannot be shrunk to a point because it encircles 
a hole in the superconductor, we write v =ae’S, 
where a and S are real; then In(¥/¥")=2iS. Since 
W must be single valued, S can change only by 
2mn where n is an integer, so Eq. (5) yields 


&, =n(nhc/e). (6) 


This is the result we sought. It is a characteris- 
tic feature of superconductors which distinguishes 
them from normal conductors even when they are 
pervaded by magnetic fields, as is the case for 
specimens which are thin compared to the pene- 
tration depth :. From (6) it follows that the flux- 
oid remains constant in time in the superconduct- 
ing state, since it must change continuously and 
n is an integer. When J=0 along C, Eq. (2) shows 
that $c is just the magnetic flux through C. Then 
the flux itself is a multiple of the natural unit. 
This is the fact which was verified experimentally. 
Bardeen’s result can be deduced easily from Eq. 
(6) by considering a thin cylindrical tube of radius 
y and thickness d$)\. We assume that the current 
J is constant throughout the tube and in the circum- 
ferential direction and that the density ||? is 
also constant. Then the total current per unit 
length of tube is d|J| and it produces a flux $= 
4n*y7d|J|/c through the hole. The line integral 
of A in Eq. (2) is just # if C is a circle of radius 
y concentric with the axis, so Eqs. (2) and (6) 
yield 


+ 2nvcA |J| =n(ahc/e). (7) 


Upon introducing the penetration depth \ = 
(mc? /8ne? |¥|?)? and eliminating |J| in terms 
of @, we obtain from Eq. (7) 


oan (1 -®) 4 (8) 
e rd 

This is just Bardeen’s result. Our derivation 
shows that an experimental verification of Eq. (8) 
would verify the quantization of the fluxoid in nat- 
ural units. The penetration depth , becomes in- 
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finite as the temperature T approaches the criti- 
cal temperature T, from below. Therefore when 
T is sufficiently near T, and d is finite, » is much 
greater than d and so Eq. (8) applies. This shows 
that the unit of flux quantization tends to zero as 
T tends to T,. Since |J| is proportional to @ it 
also tends to zero and so does the density | ¥/?. 
However, the ratio |J|/|¥|?, which occurs in the 
fluxoid, tends to nefii/mr. 

When the tube is not thin compared to A, we have 
obtained the quantum of flux through the hole as 
well as the magnetic field as a function of radial 
distance from the axis by solving Maxwell’s equa- 
tions for A with J given by Eq. (1). By using these 
results we have computed the Ginzburg-Landau 
expression for the free energy of the equilibrium 
states. It is the sum of a surface energy propor- 
tional to n? and the free energy in the absence of 
a magnetic field, which is independent of the quan- 
tum number m. When 7 is small and the tube is 
thick compared to \, the surface energy is negli- 
gible. In any case these equilibrium values must 
be interpreted as local minima with respect to 
small deviations from equilibrium. 


We wish to thank Professor C. N. Yang for a 
stimulating discussion of his microscopic theory. 
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ELECTRICALLY INDUCED SHIFT OF THE F"® RESONANCE FREQUENCY IN MnF,1 


P. S. Pershan and N. Bloembergen 
Gordon McKay Laboratory, Harvard University, Cambridge, Masschusetts 


(Received August 2, 1961) 


The predicted linear effect of an applied elec- 
tric field on the magnetic hyperfine interaction’ 
has been found experimentally in antiferromag- 
netic MnF,. We believe this is the first such ob- 
servation, although the quadratic variation with 
E of the hyperfine splitting in the Cs atom has 
been reported earlier.” 

The F*® resonance,® which occurs at 159.970 
Mc/sec without external magnetic field at 4.2°K, 
was detected with the spectrometer built by Ku- 
shida for the investigation of the pressure depend- 
ence of the same resonance.‘ 

A plane parailel slab, 1 mmx5mmx20mm, was 
cut from a single crystal of MnF, with the [100] 
direction perpendicular to the broad face. Metal 
electrodes were applied to the broad face by sil- 
ver paint that was fired at 1000°F. Contacts were 
then attached with indium solder. A dc voltage up 
to 9700 volts could be applied from a battery of 
dry cells. It was necessary to use rf chokes to 


isolate the high-voltage line from the tank circuit. 
At the highest field, E = 0.97 10° volts/cm, the 

F’® resonance is completely resolved into two com- 
ponents. In Fig. 1 the recordings of the absorption 
derivative are reproduced for zero field and three 
values of the applied electric field. They are com- 
patible with a linear absolute shift Av; =5 kc/sec 
for each of the four F’® nuclei in the unit cell, if 

E =10* volts/cm is applied in the [100] direction, 

as shown in Fig. 2. The separation of the two com- 
ponents is of course twice the shift. 

For an arbitrary direction of E the line should 
split into four components. Only the component of 
E pointing from each F"* nucleus to the nearest 
Mn** ion is effective in shifting the resonance of 
that nucleus. If E =10* volts/cm is applied along 
a [110] direction, half the F** nuclei will keep an 
undisplaced resonance, while one quarter will be 
displaced by 7.1 kc/sec and the remaining quarter 
by -7.1 kc/sec. 
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FIG. 1. The absorption derivation of the F"® reso- 
nance in MnF, at 4.2°K in zero magnetic field. An 
electric field of variable strength is applied along the 
[100] crystallographic axis. 











This experimental result is an order of magni- 
tude smaller than the admittedly crude theoreti- 
cal estimate,’ which was based on the assumption 
of an isotropic binding force of the F~ ion in the 
(001) plane. Inspection of the nuclear distances 
in MnF, and the ionic radii of F~ and Mn** re- 
veals that the F~ ion at the position (u,u, 0) is 
most tightly bound with the Mn** ion at the origin, 
and somewhat weaker with the Mn** ions at 
($,4,+4). Its interaction with F~ ions and farther 
Mn** ions should be much weaker. The binding 
is therefore highly anisotropic. When one uses 
the known force constants of the MnF and MnF, 
molecules,® the change in the nearest Mn- F dis- 
tance Avy, could be nearly an order of magnitude 
less than estimated earlier. The assumption that 
u does not change with hydrostatic pressure may 
also have to be revised. If the internuclear dis- 
tance ry is most difficult to change, u would in- 
crease somewhat under compression. This would 
make the values for (8A/8r); yy quoted by Benedek 
and Kushida‘* too large. A revision of the covalent 
effects, utilizing the refinements pointed out by 
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FIG. 2. The absolute shift of the F'® resonance versus 
applied electric field strength along [100]. The “elec- 
tric doublet splitting” is twice the shift and amounts to 
1 eps/(volt/cm). 





Freeman and Watson, and a determination of the 
internal electric field at the bond position, is also 
necessary to obtain agreement with experimental 
observation. 

Finally it is pointed out that the observed effect 
may be used to study the presence or absence of 
antiferromagnetic domain walls. Néel introduced 
this concept, and several interesting theoretical 
discussions have appeared,®»® but experimental 
evidence is difficult to obtain. Consider a crystal 
of MnF, with an external magnetic field H, (>25 
gauss) along the c axis, [001]. Two F'® magnetic 
resonance lines will be observed. Apply an elec- 
tric field along the [110] direction. If the crystal 
is a single domain, only one component of the 
“magnetic doublet” will be split further. Work 
on this problem is in progress. 
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COMPARISON OF THE BETA SPECTRA OF B” AND NT 


T. Mayer-Kuckuk* and F. C. Michel 


California Institute of Technology, Pasadena, California 
(Received August 3, 1961) 


The observed near equality of the vector beta 
coupling constant given by the decay of O” and 
the muon decay coupling constant led Feynman 
and Gell-Mann! to postulate the conservation of 
the vector part of the beta-decay current. Gell- 
Mann’? pointed out that a consequence of this con- 
served vector current (CVC) theory is an addi- 
tional coupling in the beta interaction analogous 
to the anomalous magnetic moment of the nucleon 
(“weak magnetism”). The interference between 
this term and the Gamow-Teller transition ampli- 
tude for a 1++0+, AT =1 beta transition then re- 
sults in a deviation from the allowed shape of the 
B™ and N” spectra. The magnitude of the effect 
can be calculated?’* from the rate of the analo- 
gous M1 gamma transition in C’*. If S(E) is the 
spectrum of each transition, divided by the cor- 
responding allowed Fermi spectrum, then the 
CVC theory predicts 


S(E, B**) /S(E, N’*) = const[1+ (A + 6A)E]f(E), 


where the constant A arises from this interfer- 
ence term. The quantity f(E) is the inner brems- 
strahlung correction. The expected value for A 
is (1.33+ 0.15) % per Mev with an electromagnetic 
correction* 6A of (-0.25+0.15)% per Mev. The 
uncertainty in A comes from the uncertainty in 
the C’* gamma-decay rate and the uncertainty 
in 6A is an estimate of the error from using the 
shell model. Without the CVC theory hypothesis, 
A is estimated to be roughly 5 times smaller.*® 
The beta spectra of B” (E,,,,, = 13.369 Mev) 
and N?? (E max = 16.43 Mev) were analyzed with 
the iron-free single-lens magnetic spectrometer 
described by Hornyak et al.®° The baffle system 
was modified for the high electron energies in- 
volved. B* (20.3 msec) was produced by bom- 
barding targets of natural boron with 1.65-Mev 
deuterons, and N*? (12.5 msec) was produced 
by bombarding targets of enriched (96%) B*° 
with 2.75-Mev He® ions. Both beams were pulsed 
at 60 cps at the ion source of the electrostatic ac- 


celerator. The targets consisted of approximately 
0.3 mg/cm’ of boron deposited by cracking dibo- 
rane on foils of thickness 0.5 mg/cm? and 3 mg/ 
cm? for B* and N”, respectively. The beta rays 
were detected in a 10-mm thick anthracene scin- 
tillator. This is sufficiently thin (approximately 
2.3-Mev range for electrons) to give nearly iden- 
tical pulse-height spectra for all energies in the 
range studied (5-13 Mev). The use of a thin scin- 
tillator allows amplifier gain and discriminator 
settings to be left constant, reduces background 
as a result of its small volume, and produces a 
negligible small pulse tail extending into the noise 
background.” The activity produced on the target 
was monitored by the reaction protons from B*%(d, 
p)B™ for B” and B*°(He®, p)C for N’. A silicon 
p-n junction counter mounted in the spectrom- 
eter near the target served as a proton detector. 

The spectra were sampled at 0.5-Mev intervals 
in the range 5 to 10.5 Mev for B” and 5 to 13 Mev 
for N*”, the data being collected at alternating 
high- and low-energy points to minimize effects 
of target deterioration and instrumental drifts. 
The lower limit of 5 Mev was dictated by uncer- 
tainties in the branching ratios to excited states 
of C’*; the upper limit was maintained well below 
the spectrum end points to minimize uncertainties 
arising from calibration errors and background 
subtraction. The background was studied at ener- 
gies above the end points of the beta spectra, at 
zero field, and by using target backings without 
the boron. The neutron background was separate- 
ly investigated with a beryllium target. Typical 
background corrections at 8 Mev amounted to 
about 1.7% for B*” and 4.2% for N’*. A B™ spec- 
trum taken with the relatively thicker backing 
used for the N** showed no deviation from the 
other B’ spectra. 

The observed spectra are corrected for branch- 
ing®’® to the 4.4- and 7.6-Mev excited states of C” 
and for inner bremsstrahlung.* The N” branch 
to the 7.6-Mev state is assumed to be 3.5%. The 
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small inner bremsstrahlung correction has been per Mev and the N” energy dependence is (0.17 
applied entirely to S(E, B™). The Fermi plots, + 0.20) % per Mev giving (1.13 + 0.25) % per Mev 
not corrected for branching, are given in Fig. 1. for A +6A. 

The ‘B*? and N” matrix elements, normalized to This result appears to be in agreement with the 
1.00 at 8 Mev, are displayed in Fig. 2. The B™” CVC-theory prediction of (1.08+ 0.22) % per Mev. 
plot represents one typical run, and the N” plot An important uncertainty is the 7.6-Mev branch- 


is the data of 8 runs. We find the B’” matrix ele- ing ratio in the N™ decay. If, for example, we 
ment has an energy dependence of (1.304 0.10)% have overestimated this branching ratio by 10%, 





FIG. 1. Fermi plots 
for B'* and N’*, The ob- 
served end-point energies, 
indicated in the figure, 
are in agreement with the 
values from reaction data 
[F. Ajzenberg-Selove and 
T. Lauritsen, in Landolt- 
Bérnstein Tables (to be 
published)] of 13.369 Mev 
for B!* and 16.43 for N"’, 
The curvature of the N” 
plot due to branching can 
easily be seen. In the B” 
plot this effect is overcom- 
pensated by the CVC shape 
factor. Branching end 
points are indicated by 
arrows. 
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FIG. 2. Plots of the shape factors S(E,B"*)/f(E) and 
S(E,N"?) normalized to unity at 8 Mev. Statistical errors 
and least-squares fits are given. Branching end points 
are indicated by arrows. Slopes are 1.29% per Mev 
for B'* and 0.17% per Mev for N!?, 


our value for A + 6A should be increased by 0.12% 
per Mev. It is interesting that the relatively 
strong energy dependence of the B*” matrix ele- 
ment accounts for the failure of Hornyak et al.’° 
to observe a B™ branch to the 4.4-Mev state of 
Cc, and in fact gives a Fermi plot that is slightly 
concave downwards. 

Previously, Nordberg, Morinigo, and Barnes™ 
compared the asymmetries of the beta-alpha an- 
gular correlation in the decays of Li® and B® and 
also found agreement with the CVC theory. A 
small discrepancy remains, however, between 
the muon lifetime and the Fermi coupling con- 
stant derived from the O“ decay.” Possibly this 
discrepancy is due to a charge dependence of the 
internucleon potential.** 


The authors are indebted to Professor T. Lau- 
ritsen, Professor C. A. Barnes, Professor W. A. 
Fowler, Professor H. A. Weidenmiiller, and 
Professor F. Boehm for valuable discussions. 
Two earlier experiments in this laboratory by 
H. Hilton and V. Soergel, and H. Hilton and C. 
Van der Leun provided much of the groundwork 
which made the present experiment possible. 

One of us (Mayer-Kuckuk) was supported by 
Deutsches Bundesministerium ftir Atomkernener- 
gie. 
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DOUBLE GAMMA EMISSION IN THE 6.06-Mev MONOPOLE TRANSITION OF O**T 


S. Gorodetzky, G. Sutter, R. Armbruster, P. Chevallier, P. Mennrath, F. Scheibling, and J. Yoccoz 
Institut de Recherches Nucleaires, Strasbourg, France 
(Received July 27, 1961) 


In a continuation of our work on double gamma- 
ray emission in electric monopole transitions we 
have undertaken a study of the 6.06-Mev transi- 
tion in O**. This case was considered theoreti- 
cally by Oppenheimer and Schwinger’ in 1939, 
and then by Dalitz? and by Cameron.’® Very re- 
cently, Yoccoz* has proposed a more complete 
theory of this problem and gives as preliminary 
result for the ratio of the probabilities of double 
gamma-ray emission and pair creation: 3x107% 
<Iy,/T, <8x10™, At present, more detailed 
calculations are being undertaken. 

The 6.06-Mev level was excited by the reaction 
F'°(p, a)O**, at the 1880-kev resonance. An exci- 
tation curve® for nuclear pair emission had shown 
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that this resonance was the most advantageous one 
for a study of the monopole transition. 

The method of summing the energies of the two 
gamma radiations, previously used in an attempt 
to detect double gamma-ray emission from the 
1.75-Mev transition in Zr™,® was not suitable for 
the 6.06-Mev transition in O'*, Indeed for the 
gamma-ray energies occuring in O** the total ab- 
sorption rate in the crystals used is much lower 
than in Zr®, We have thus been led to the use of 
(a-y-y) triple coincidence (27=5x10~*sec), 
gating a multichannel analyzer which records the 
energy spectrum of the alpha particles (Fig. 1). 
The gamma-rays were detected by two 4-in. x4-in. 
Nal crystals at 180°, and the alpha particles by a 


Ao 





FIG, 1. (a) Ep = 1880 kev, 
eo 712 direct alpha-particle spectrum. 


«6.92 (b) Ep = 1880 kev, alpha-parti- 
nen cle spectrum gated by alpha- 
ae gamma-gamma triple coinci- 
dences. Gamma-ray window 
at 0.511 Mev (annihilation ra- 
2 im diation). Accidentals subtract- 








692 ed. (c) Ep = 1880 kev, alpha- 

2 = particle spectrum gated by 
alpha-gamma double coinciden- 
ces (E,>2 Mev). Accidentals 
subtracted. (d) Ep =1375 kev, 
alpha-particle spectrum gated 
by alpha-gamma double coinci- 
Xo dences (Ey>2 Mev). Acciden- 
tals subtracted. 
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FIG, 2. Ey =1880 kev, alpha-particle spectrum gated 
by alpha-gamma-gamma triple coincidences. Gamma- 
ray windows from 1.5 to 4.5 Mev. Accidentals not sub- 
tracted (see Table J). 


20-mm? RCA junction detector at 150° to the direc- 
tion of the incident beam. 

A preliminary experiment’ had given an upper 
limit: Pyy/T a <6x10%, 

Since this first experiment the apparatus has 
been improved by the use of a lead target cham- 
ber forming a conical diaphragm for the Nal crys- 
tals. Only gamma rays between 1.5 and 4.5 Mev 
were accepted by the circuit. 





Ny,N Cc 


At E,=1880 kev a series of 166 measurements, 
each of 30 minutes duration, were made. Figure 
2 shows the spectrum obtained. We have verified 
that the distribution of the coincidences in time 
follows the Poisson law. 

Since the measurement of the accidental coinci- 
dence rate would be excessively time consuming 
because there are several possible kinds of such 
coincidences, we calculated® the rate from the 
various counting rates and the coincidence resolv- 
ing time. Measurements made at the 1375-kev 
resonance [predominantly producing the 6.13-Mev 
level—see Fig. 1(d)] seem to show that the method 
is correct. The probability of double gamma-ray 
emission from this level is theoretically very 
small. Under the same experimental conditions 
as at 1880 kev (in particular, the same counting 
rate) no real events were observed; all the coinci- 
dences obtained could be attributed to accidentals. 
Furthermore, this method has been verified at the 
1880-kev resonance by measuring the accidental 
coincidences of different kinds that occur with a 
greater beam intensity. We estimate that the 
maximum error in the calculation of accidentals 
is 20%. 

Results obtained at the two resonances are 
shown in Table I. 

The experimental results are related to the 
double gamma-ray emission probability by the 
following formula: 


7 
Eo /140 
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where E,= 6.06 Mev,(Q2n)gy= the mean value of the 
product of the solid angle and the efficiency of a 
crystal, e€= the efficiency of the coincidence cir- 
cuit (here «=1), C=0.89=a factor due to the 

(1+ cos?@) correlation between the two El gamma 
rays detected, and a=0.87= the mean transmis- 
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sion of gamma rays by the absorber which stops 
the electrons from nuclear pair creation. 

The product (2n)gy was measured directly at 
the 1375-kev resonance from the ratio 


N coine'*~”) 6.13" 06.13” 


Table I. Alpha-gamma-gamma measured coincidences and calculated accidentals. 
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*% 6.92% 7.12 “6.13 a 
1375 kev Measured coincidences 40 55 
Calculated accidentals 43 52 
1880 kev Measured coincidences 146 55 
Calculated accidentals 147 15 
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taking into account the a-y angular correlation. 
We obtained (Qn), y€=3.66x10™. This result 
was verified by a measurement made with a cal- 
ibration source of Co” corrected to our energy 
by the use of curves given by Wolicki et al.° 

Using the energy spectrum of double gamma- 
ray emission, which has the form E*(E,-£)*,!°"'*»*, 
we have calculated the integral in the denominator 
of Eq. (1) and find 


4.5 Mev 
[Eo(E)P [€)-£) PE, -£) Pak = 0.24, 


(2) 


where p(E) is the form factor of the energy spec- 
trum of a gamma ray of energy E in a 4 in. x4 in. 
Nal crystal. The term 140/E,” is the normaliza- 
tion factor. The number No is measured direct- 
ly during the experiment. In a separate experi- 
ment the ratio Ngo/Nqm has been measured as 
2.8. The factor 1/2 is due to the indistinguishibil- 
ity of the two gamma rays. 

Since approximate calculations have shown that 
the internal and external bremsstrahlung coming 
from the pairs are negligible,’* we believe that 
the true coincidences observed can be attributed 
to double gamma-ray emission. Insofar as no un- 
known phenomena of the same order of magnitude 
are taking place, our experimental result leads 
to the value 


140 


7 
Fo Ji.5 Mev 


r /P_=(2.5+1.1)x10%, 
vy 7 


In addition, the measurement made at 1375 kev 
allows us to give an upper limit for the double 
gamma-ray decay of the 6.13-Mev £3 transition 








of O#*: 
r /r <5x104, 
» ie 
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MAGNETIC MOMENT OF N*T 


M. Posner, J. L. Snider, A. M. Bernstein, and D. R. Hamilton 
Palmer Physical Laboratory, Princeton University, Princeton, New Jersey 
(Received August 15, 1961) 


Uncertainties in the form of the magnetic dipole 
moment operator and the exact structure of nucle- 
ar wave functions make the theoretical interpreta- 
tion of nuclear magnetic moments difficult.’ In 
addition to the usual orbital and spin contributions 
to the moment one must consider the effects of 
relativity’ and of the currents due to the exchange 
of charged 7 mesons between nucleons (exchange 
moments).* As a consequence of the assumption 
of the charge independence of nuclear forces the 
wave functions of mirror nuclei are the same. A 
comparison of the magnetic moments of a mirror 
pair could, in principle, yield information about 
the structure of the operator. In particular, Sachs 
has pointed out that the sum of the moments of a 
pair of mirror nuclei should be independent of 
mesonic current contributions.*® 

With these ideas in mind we have measured the 
spin and magnetic moment of 10-minute N** by 
the atomic beam magnetic resonance method. Be- 
cause of the short half-life a continuous-flow 
source was used. The radioactivity was produced 
by the bombardment of N**, in melamine powder, 
with 18-Mev protons. The resulting N*’, in gas- 
eous molecular form, was flushed with neon from 
the cyclotron target to the atomic beam machine. 
There a microwave discharge dissociated the 
molecules to produce a beam of neutral N** atoms. 
These atoms were then passed through a machine 
of a new design, consisting of a six-pole focusing 
A magnet,* a C magnet,® and a two-pole, constant- 
gradient, B magnet. The atoms were collected by 
means of a chemical reaction on titanium foils 
maintained at a temperature of 1050°C. After ex- 
posure to the beam these foils were removed from 
the apparatus and the positron activity was counted 
in standard scintillation counters. 

The optics of the machine allowed us to place 
two collector foils at the B-magnet exit, one of 
which monitored the beam while the other collected 
primarily atoms which had undergone resonance 
transitions. The ratio of the count rates on these 
foils is low off resonance, increases at resonance, 
and is independent of the size of the beam. This 
technique enabled us to observe resonances with 
beams which were not completely steady. 

The N** atoms were produced mainly jn the ‘S,, 
atomic ground state, but resonances were also 
seen in the metastable *D,, and *P,, states. A 


nuclear spin of $ was measured by observing 
resonances in the Zeeman region (where the reso- 
nance frequencies vary linearly with magnetic 
field) in the atomic states mentioned above.® Sev- 
eral resonances were observed at each of a num- 
ber of magnetic fields, these being consistent 
only with a nuclear spin of 5. 

To measure the value of the hyperfine structure 
separation, Av, the *S,., F =2, mp =1 to F =2, 
my =0 transition (assuming » <0) was followed 
out of the Zeeman region to a field of 22 gauss. 
The shifts of the resonance frequency from line- 
arity with field give 


Av =33.2+0.3 Mc/sec. 


Combining this with the value of Av and py for N**,” 
we obtain 


| (N#%) | =0.321+ 0.003 nm. 


One would assume that the sign of the moment is 
negative as in N*®. 

In Table I we list the data now available on the 
magnetic moments of pairs of mirror nuclei. 
(This table does not include the neutron and pro- 
ton whose moments are basic ingredients of the 
single-particle or Schmidt moments.) For both 
pairs of nuclei the moments lie “outside” the 
Schmidt moments by amounts which are approx- 
imately equal and opposite. From a careful con- 
sideration of the sum of the moments of the mass- 
three nuclei, Sachs has shown’ that the moments 


Table I. Magnetic moments of mirror pairs of nuclei. 
KE is-the experimental moment, us the Schmidt single- 
particle moment. The values of ys for He® and H® are 
the same as the free neutron and proton moments, upon 
which the Schmidt moments are based. All values of BE 
except that for N"’ are taken from Strominger et al. * 








Nucleus m1 E He (u ne m 5) 
He® -2.127 -1.913 -0. 214 
H® 2.979 2.793 0. 186 
Sums 0.852 0.880 -0. 028 
c3 0.702 0.638 0. 064 
nis (-)0.321 -0. 264 -0. 057 
Sums 0.381 0.374 0.007 





3D. Strominger, J. M. Hollander, and G. T. Seaborg, 
Revs. Modern Phys. 30, 585 (1958). 
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cannot be explained by any reasonable mixing of 
LS states and that therefore these results must 
be due to mesonic current contributions of 0.27 
nm. 

For the A =13 mirror pair, calculations per- 
formed using the intermediate-coupling shell 
model and the usual magnetic moment operator 
show that if there is a sufficient deviation from 
jj coupling toward LS coupling, the magnetic mo- 
ments can lie outside the Schmidt limits.*"° An- 
alysis of the available experimental data, not in- 
cluding the magnetic moments or M1 transition 
rates, indicates that the coupling is intermediate 
between LS and jj, but is not of sufficient preci- 
sion to determine the exact coupling.*° Therefore, 
the prediction of the model for the magnetic mo- 
ments cannot be precisely determined at present. 

It should be pointed out that the Hamiltonian 
used in these calculations includes exchange and 
velocity-dependent (spin-orbit) interactions which 
change the form of the current operator. Since 
the proper form of this operator is not known, 
these terms have not been taken into account. 
Therefore, the calculation of the magnetic mo- 
ments and magnetic transition rates is not self- 
consistent." 

The sum of the moments of N** and C"* as cal- 
culated in intermediate coupling is independent 
of the coupling and is equal to the sum of the 
Schmidt moments to within 0.04 nm,” in agree- 
ment with the experimental results. From this 


agreement one can conclude that the mesonic cur-- 


rent contributions in N** and C*’, if they exist, 
obey the mirror principle (i.e., they are equal 
and opposite in a pair of mirror nuclei).*’** 


A more detailed discussion of this work will 
be published subsequently. 





+ 
This work was supported by the U. S. Atomic Energy 


Commission and the Higgins Scientific Trust Fund. 
‘For a review of the subject see R. J. Blin-Stoyle, 


Theories of Nuclear Moments (Oxford University Press, 





New York, 1957). 

21. Rosenfeld, Nuclear Forces (North-Holland Pub- 
lishing Company, Amsterdam, 1949), p. 417. 

3R. G. Sachs, Nuclear Theory (Addison-Wesley Pub- 
lishing Company, Inc., Reading, Massachusetts, 1953), 

‘R. L. Christensen and D. R. Hamilton, Rev. Sci. 
Instr. 30, 356 (1959). 

5L. S. Goodman, Rev. Sci. Instr. 31, 1351 (1960). 

63. L. Snider, M. Posner, A. M. Berstein, and D. 
R. Hamilton, Bull. Am. Phys. Soc. 6, 224 (1961). 

"L. W. Anderson, F. M. Pipkin, and J. C. Baird, 
Jr., Phys. Rev. 116, 87 (1959). 

8For an introduction to the intermediate-coupling 
shell model see J. P. Elliot and A. M. Lane, Handbuch 
der Physik , edited by S. Fliigge (Springer-Verlag, 
Berlin, 1957), Vol. 39, p. 241. 

*D. Kurath, Phys. Rev. 101, 216 (1956). 

104. M. Lane, Proc. Phys. Soc. A68, 197 (1955), 
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‘In principle, it seems reasonable that these currents 
should be equivalent to the mesonic exchange currents 
since the nuclear force is believed to derive from the 
interchange of mesons. 

2D. Kurath, Phys Rev. (to be published). We thank 
Dr. Kurath for sending us a preprint in advance of pub- 
lication. 

13We have not discussed the relativity contribution to 
the magnetic moment as it is expected to be small com- 
pared to the deviation from the Schmidt moment. See 
reference 2. 











IDENTIFICATION OF DOUBLET STATES AT 5.16 Mev in B*° 


E. L. Sprenkel, J. W. Olness, and R. E. Segel 
Aeronautical Research Laboratory, Dayton, Ohio 
(Received August 15, 1961) 


The isotopic spin analog of the first excited 
states in Be*® and C*® should be located at an ex- 
citation energy of about 5.1 Mev’ in the self-con- 
jugate nucleus B’®. In this region three levels 
have been reported previously at 4.77, 5.11, and 
5.16 Mev, of which the first two levels may readi- 
ly be eliminated as candidates for the J” = 2", 
T=1 isotopic spin analog. The energy displace- 
ment of the 4.77-Mev state, from the expected 
position at 5.1 Mev, would be exceptionally large 
and the broad width (['¢_m, =1.2 kev) of the 5.11- 
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Mev state makes an assignment of T=1 very un- 
likely for that state. Recently Meads et al.* have 
shown, from a comparison of proton and deuteron 
(about 10-Mev incident energy) inelastic scatter- 
ing by B*°, that the isotopic spin of both the 4.77- 
and 5.11-Mev states is T=0. 

The state at 5.16 Mev has presented a puzzle. 
Investigation of the 5.16-Mev level when fed by 
2.40-Mev gamma-ray transitions from the 7.56- 
Mev level in B*° yields results which are not sup- 
ported when the 5.16-Mev level is produced by 
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other means. From a study of the angular distri- 
butions of the gamma rays emitted from the 5.16- 
Mev level, produced in the Li®+a@ reaction, 
Meyer -Schiitzmeister and Hanna‘ found that the 
possible J” assignments were 1* or 2*, the decay 
scheme favoring the latter assignment. The ap- 
parent absence of transitions to the J” = 0°, T=1 
state at 1.74 Mev would appear to rule out a J” 
=1*, T=0 assignment to the 5.16-Mev level on 
the basis of Morpurgo’s rule.® Since aJ"=1", 
T=1 assignment cannot be reconciled with the 
T=1 system of levels in the mirror nuclei Be’® 
and C!°, the proper assignment for the 5.16-Mev 
level would be expected to be J"=2*, T=1. Re- 
cently the 7=1 assignment has been supported by 
Meads et al.* as a result of their comparison of 
proton and deuteron inelastic scattering. On the 
other hand, the strong 2.40-Mev gamma ray emit- 
ted by the J7=0* state at 7.56 Mev in B’° indicates 
that the assignment of the 5.16-Mev level® popu- 
lated by this transition must be J=1. 

The intrinsic c.m. width of the 5.16-Mev level 
of B*° formed in the Li®+a reaction has been re- 
ported by Hanna’ to be less than 500 ev. Ina pre- 
liminary investigation of Be®(d,n)B*° reactions, 
Chase and Warburton® report both neutron and 
gamma-ray thresholds for the 5.16-Mev level 
with r,>0.1f. For Be*(p,y)B*° reactions, Mey- 
erhof, Tanner, and Hudson® have shown, from 
the spectra in coincidence with capture gamma 
rays to the 5.16-Mev level, that r,,~0.1f for 
this state. Similar results were obtained by 
Meyerhof and Chase’® for the gamma rays emit- 
ted from the interaction of 2.8-Mev deuterons 
with a thick beryllium target. On the other hand, 
Sprenkel and Daughtry” did not observe gamma 
rays in coincidence with 2.40-Mev transitions 
from the 7.56-Mev level and placed a limit of 
r,<0.01P for the 5.16-Mev state. In order to 
explain the disagreement of their results with 
those presented above, Sprenkel and Daughtry” 
postulated the existence of doublet states at 5.16 
Mev in B®, 

Since the doublet state that is fed by 2.40-Mev 
gamma rays is expected to have aJ=1, T=0 as- 
signment, the decay of this state should proceed 
mainly through the transition to the J7=0*, T=1 
State at 1.74 Mev in B'*°. The expected predomi- 
nance of this transition (3.42 Mev) arises from 
the selection rules for the strong reduction in 
Strength of E1 and M1 transitions between levels 
of the same isotopic spin in self-conjugate nuc- 
lei.®»42 Accordingly a search was made for a 3.42- 
Mev gamma-ray transition that should result from 


the decay of a J=1 doublet member produced by 
bombarding Li® with alpha particles of about 1.17 
Mev. Since the known member of the 5.16-Mev 
doublet decays by 3.01-Mev (64%), 4.44-Mev 
(29%), and 5.16-Mev (7%) transitions,* the detec- 
tion of a peak corresponding to the full energy of 
a 3.42-Mev gamma ray is rendered difficult, be- 
cause the peak lies in the same spectral region as 
the second escape peak of the 4.44-Mev transition. 
In order to increase the sensitivity for the detec- 
tion of 3.42-Mev gamma rays, spectra were taken 
in coincidence with either of the other members 

of the 3.42-1.02-0.72 Mev cascade. Two Nal(T1) 
crystals (8 in. diameter x8 in. long and 5 in. di- 
ameter x5 in. long) were placed at 90° with re- 
spect to the alpha beam and at a distance of about 
3/4 in. from the Li® target. To minimize escape- 
peak formation, the larger Nal(T1) crystal was 
used for the detection of the coincidence spectrum. 
The fixed channel, for the smaller crystal, ex- 
tended from 0.6 to 1.2 Mev for selection of the 
0.72- and 1.02-Mev cascade gamma rays. Coinci- 
dence spectra were obtained at alpha-energy inter- 




















Li*(a, y)B° 

100+ E,* 1.01 Mev 
g °F 
cs Lge ; Pe 
= £ Sheik Lani oa Wi Sesh en edeanatennd 
E 
Ss , 
oe 100} “ E, = 1.31 Mev 
o . es 3 
as e bag ~ ‘ ’e. i 
3 Tenet: fe . 
3 iii ° nngat san, 
w 4 > iene tans 
o 
D wok, E, = 1.62 Mev 
z 7 
sa i”. 
Oo sol- " 

RAIN as ¥ 
Seth Se. oot 
SP Sapa 
Ey (Mev) 


FIG, 1. Gamma-ray coincidence spectra in the re- 
gion of the 1.20-Mev resonance in Li*(a,73,4,)B'®. The 
middle spectrum also shows gamma rays from the pre- 
viously known 1.175 - Mev resonance. 
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vals of 75 kev using an enriched Li® target with an =1.31 Mev (Fig. 1). For positive identification 
average thickness of 50 kev. A peak correspond- of the 3.42-1.02-0.72 Mev cascade chain, a 

ing to a 3.42-Mev transition was seen to appear at spectrum was taken in coincidence with the “3.42- 
about E ,=1.0 Mev, increase sharply in intensity Mev” gamma ray detected in the large crystal. 

to reach a maximum at about E a" 1.2 Mev, and The resultant spectrum showed the presence of 
then gradually decrease in intensity (Fig. 1 and only two gamma rays at 0.72 and 1.02 Mev, which 
Fig. 2). The position of the peak shifted ina reg- provided identification of the expected cascade 
ular manner (Fig. 1), showing a change ingamma-_ chain. 





ray energy from 3.36 Mev at E,, =1.0 Mev to 3.70 The “3.42-Mev” gamma-ray yield curve (Fig. 
Mev at E,,=1.6 Mev. Peaks corresponding to 2) was fitted with the single-level expression: 
3.01- and 4.44-Mev transitions from the previ- 
ously known member of the doublet were found r at » 

= 2 ’ 
to be very intense in the coincidence spectrum "= ‘Ee +a EF *T, 7a? (1) 


taken at E= 1.19 Mev. Because of the nonuni- 
form thickness of the target, the 3.01-Mev peak 
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parameters is given by Sachs.'* Using an inter- 
action radius of 4.94x107'* cm and assuming s- 
wave formation, the theoretical curve satisfactor- 
ily fit the experimental data when E,(lab)=E) + A,° 
=1,210+ 0.035 Mev (B*° =5.18 Mev) and (lab) 
=0.340+ 0.050 Mev. The ~200-kev (c.m.) breadth 
of the state agrees with the report that the 2.40- 
Mev gamma ray feeding the state appears to have 
an observable energy spread."'»** The width cor- 
responds to about 86% of the single-particle limit 
for s-wave (a) formation and 190% for p-wave for- 
mation. Be*%(d,n)B*° angular distribution data’® 
imply a negative-parity state in the 5.1-Mev re- 
gion for which the newly found state offers the only 
possibility. However, the large width measured 
here favors s-wave formation although there is 
enough uncertainty in the resonance parameters 

to allow p-wave formation of approximately single- 
particle width, 

Comparison of the thick-target yields for the 
3.01- and 3.42-Mev gamma rays, using the previ- 
ously determined* 3.01-Mev gamma-ray yield, 
provided the value [3 42=0.06+ 0.03 ev. This 
radiation width corresponds to |M|?~0.06 which 
is well within the predicted’® range of an allowed 
M1 transition. An upper limit of ~1/3 the intensi- 
ty of the 3.42-Mev transition was placed on any 
other gamma ray emanating from this broad 5.18- 
Mev state. The branching ratio f/f ~3x10~-" 
is consistent with the upper limit of 10 estab- 
lished by Sprenkel and Daughtry." 

The identification of another state in the 5.1- 
Mev region of B’° finally allows a solution of the 
puzzle regarding the isotopic spin analog of the 
first excited states in Be’® and C’®, The 7.56- 


Mev level decays by 2.40-Mev gamma rays to 
feed the J=1, T=0 broad level centered at 5.18 
Mev, leaving the previously known narrow 5.16- 
Mev level for the long sought J7=2*, T=1 state. 
These assignments allow reconciliation of all con- 
flicting results from the previous experiments 
described above. 
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EVIDENCE FOR A T=0 THREE-PION RESONANCE* 


B. C. Magli¢, L. W. Alvarez, A. H. Rosenfeld, and M. L. Stevenson 
Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 
(Received August 14, 1961) 


The existence of a heavy neutral meson with 
T=0 and J=1~ was predicted by Nambu’ in an at- 
tempt to explain the electromagnetic form factors 
of the proton and neutron. Chew’ has pointed out 
that such a vector meson should exist on dynami- 
cal grounds as a three-pion resonance or a bound 
state. Such a particle is also expected in the vec- 
tor meson theory of Sakurai® and, as a member 
of an octet of mesons, according to the unitary 
symmetry theory‘; and for other reasons.*° We 
will refer to it as w. 

Previous searches® for w have primarily been 
confined to the mass region m,,<3u, with py =the 
pion mass, where only the following radiative de- 
cay modes are allowed: w+7°t+y, w+2m°+y, 
and w+7*+7"+y. The w cannot decay into two 
pions. 

The present search was made assuming m,, 
>3m,, where the decay 


wont+n-+7° (1) 


is possible.” We have searched for such a 3-pion 
decay mode by studying the effective mass distri- 
bution of triplets of pions in the reaction 


Prpontents+n-+07+ 7°, (2) 


We have measured 2500 four-prong events pro- 
duced by antiprotons of 1.61 Bev/c in the 72-inch 
hydrogen bubble chamber.® The c.m. energy is 
2.29 Bev. Upon fitting these 2500 four-prong 
events by using our kinematics program KICK, 
800 four-prong events had a x? <6.5 for hypothe- 
sis (2) and would not fit the hypothesis that no 7° 
was produced (610 of these 800 had a y?< 2,5). 

The 800 four-prong events must have some 
small contamination of events in which two 7°’s 
were produced, but inspection of the “missing 
mass” distribution convinces us that it is <7%. 
Other tests confirm this low contamination. For 
example, the angular distribution of the 7° is sym- 
metric within statistics, and the momentum of the 
7° resembles the momentum distribution of the 
charged pions. 

We have evaluated the 3-body effective mass, 


M,= [(E,+E,+E) - (P, +P, + P,)?}”, (3) 


for each pion triplet in Reaction (2). Each of the 
800 four-prong events yields ten such quantities 
corresponding to the following charge states: 
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1Q1=0: a*x “n° (800x4 combinations), (4) 
1Q\|=1: a*n2*n* (800x4 combinations),  (4’) 
and 


1Q|=2: m*2*n° (800x2 combinations). (4’’) 


For each value of M, as given by Eq. (4) we can 
calculate an uncertainty 5M,, by using the vari- 
ance-covariance matrix of the fitted track varia- 
bles, which is evaluated by KICK. By using these 
6M, we have formed the resolution function of M,, 
and find that it has a half-width at half-maximum, 
T'reso}/2, equal to 8.7 Mev. However, our input 
errors to KICK allow only for Coulomb scattering 
and estimated measurement accuracy, and do not 
account for optical distortion and unknown system- 
atic errors. For example, our distributions have 
the correct shape but are too wide by a scale fac- 
tor of about 2. This suggests that our average in- 
put error is too small by about /2. Hence, our 
estimate of 5M, must be increased by about v2, 
and of T'yego)/2 to 12 Mev. We chose 20-Mev his- 
togram intervals for plotting our M, distribution. 

In Fig. 1 we have plotted the M, distributions 
for the 800 Reactions (2). Distributions 1(A) and 
1(B) are for charge combinations |Q|=1 and 2, 
respectively. The solid curves are an approxima- 
tion to phase space. 

The neutral M, distribution, 1(C), shows a peak 
centered at 787 Mev that contains 93 pion triplets 
above the phase-space estimate of 98. To contrast 
the difference between the neutral M, distribution 
and that for |Q|>1, we have replotted at the bot- 
tom of Fig. 1 both the neutral distribution and $ 
the sum of the |Q|=1 and |Q|=2 distributions. 

Figure 2 shows the M, spectra with phase space 
subtracted. The absence of the peak in the |Q| >0 
distributions determines the isotopic spin of the 


resonance, 
T =0. 
w 


The x? distribution of the events in the “peak 
region” was compared with the yx? distribution of 
the events in the adjacent “control region, ” rang- 
ing from M, 2 820 to M,<900 Mev. These distri- 
butions agree with each other, which indicates 
that the events in the peak are genuine, rather 
than being caused by some unknown background 
reaction which was misinterpreted as Reaction 
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FIG. 1. Number of pion triplets versus effective 
mass (M3) of the triplets for reaction p+p— 2m* +27" 
+7, (A) is the distribution for the combination (4’), 
|@|=1; (B) is for the combination (4"), |Q|=2; and (C) 
for (4), Q=0, with 3200, 1600, and 3200 triplets, re- 
spectively. Full width of one interval is 20 Mev. In 
(D), the combined distributions (A) and (B) (shaded 
area) are contrasted with distribution (C) (heavy line). 


(2). The missing-mass distributions in the two 
regions also agree with each other, thus support- 
ing the above conclusion. 

The peak in Fig. 2(B) appears to have a half- 
width [/2<15 Mev. This is so close to our reso- 
lution, ['yego)/2, of 12 Mev that we cannot unfold 
it without further study and at present can only 
conclude that 


M ='787 Mev, 
WwW 


and 
r/2<15 Mev. (5) 


By using the uncertainty principle, we see that 
this half-width implies a mean life t>4x10™° 
sec. Our w’s are produced with a typical c.m. 
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FIG. 2. (A) M; spectrum of the pion triplets in the 
combined distributions 1(A) and 1(B), with the smooth 
curve subtracted. (B) M; spectrum of the neutral pion 
triplets in distribution 1(C), again with the smooth back- 
ground subtracted; a resonance curve is drawn through 
the peak at 787 Mev with [/2=15 Mev. The error flags 
are VN, where N is the total number of triplets per 20- 
Mev interval before subtraction of the smooth back- 
ground curve. 


momentum of 800 Mev/c, so that in a mean life 
they travel farther than 13 f. 

We now assume that the w peak is real, and 
want to estimate how many w mesons it contains. 
As shown in Fig. 1(C), 191 triplets have M, val- 
ues between 740 and 820 Mev. (We call this the 
“peak region.”) However, these 191 triplets 
come from only 170 different four-prong events 
[i.e., 21 Reactions (2) have two values of M, in 
the peak region]. We use the charged M, distri- 
bution to estimate the background in the interval 
as 98 triplets, and then calculate a production of 
83+16 w mesons out of 800 Reactions (2); i.e., 
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(10+ 2)% of Reactions (2) proceed via 
p+pomt+n +. (6) 


Among the same 800 five-pion events, we have 
searched for—and found—the T=J=1 pion-pion 
resonance (p meson).® We found that approximate- 
ly 30% of them proceed via 


pt+rp+27+p. (7) 


We have checked whether there is any correlation 
between the observed p mesons and the w mesons. 
For each triplet inside the peak region, 740<M, 
< 820 Mev, we have evaluated the effective mass 
of the remaining 7*7~ doublet, M,. The M, dis- 
tribution is consistent with a continuum, starting 
from about 300 Mev, that has (10+ 2)% of the doub- 
lets with values of M, in the region of p, which we 
took to be 750+ 50 Mev. There is no evidence that 
the w and the p are produced in association. 

Although the masses of w and p differ by only 
35 Mev, we believe that they cannot be the same 
particle, because of their different widths (the 
['/2 for p being 40 Mev), isotopic spin, and G- 
conjugation parity—which forbids 27 + 37. 

In referring to the T=0 37 resonance as w, we 
have tacitly supposed that it is in fact a vector 
state with J=1". However, the spin and parity 
must be decided by experiment. Even if we as- 
sume the spin is <2, there are left three possibil- 
ities which are listed in Table I. A T=0 state of 
three pions must be antisymmetric in all pairs; 
hence all three pions must have different charges, 
i.e., 7°n°n° is forbidden. The matrix element of 
the m+n“n° state is conveniently analyzed in terms 
of a single pion plus a di-pion. The pions of the 
di-pion are assigned momentum Pand angular 
momentum L (in the di-pion rest frame). Then 
another pair of variables, p and I, describe the 
remaining pion in the 37 rest frame. Because the 
state is antisymmetric in any pair, L must be odd; 
henceforth, we assume L=1. Then if /=0 we 
have aJ=1" (i.e., vector) matrix element, as 





listed on the bottom line of Table I. Since three 
pions are involved, there is an intrinsic parity of 
(-1)8, so that the corresponding “meson” is not 
V, but A. 

If /=1, the matrix element can be 1+ (axial) or 
0+ (scalar) corresponding, respectively, to a vec- 
tor meson (w) or a pseudoscalar (PS) meson. 

Do we have enough data to distinguish between 
totally antisymmetric A vs S vs V matrix ele- 
ments? It is convenient to make a Dalitz plot’® 
[ Fig. 3(D) for the peak region events, 3(A) for the 
control region events] that displays the threefold 
symmetry of three pions in an antisymmetric 
state. Unit area on a Dalitz plot is proportional 
to the corresponding Lorentz-invariant phase 
space, so that the density of plotted points is pro- 
portional to the square of the matrix element. It 
is easily shown that the size of the figure is pro- 
portional to T,+ T,+ 7T,;=Q=m,, - (2mq++mzqo0)." 
Because of the finite width of the peak and the 
control regions, Q varies from event to event, so 
we use normalized variables, T;/Q. The antisym- 
metry allows the plot to be folded about any medi- 
an, so that in Figs. 3(C) and 3(B) all the data have 
been concentrated into 4 of the plot area; the sta- 
tistical distribution of the events is then more evi- 
dent. 

All three competing matrix elements, being anti- 
symmetric, must vanish where any two pions 
“touch” in momentum space. If two pions touch, 
the third must have its maximum kinetic energy 
[regions (d), (f), and (b) on the plot]. The reso- 
nance region points [Figs. 3(C) and 3(D)] seem to 
show the required depopulation at points (d), (f), 
and (b). 

More evident, however, on the plot is the fact 
that near p=0 [points (a), (c), and (e)] the densi- 
ty of peak-region points is only one half of that on 
the control plot. This is all the more suggestive 
when it is remembered that even the peak-region 
data contain only (43+ '7)% resonance events. This 
depopulation at p=0 suggests an angular momen- 


Table I, Possible three-pion resonances with T=0, J<=1. 








“Meson” Matrix element 
Type, I i, L Type, I Simple example Vanishes at: 
V,1- 1,1 A, 1+ E (Py * P,) + Eq(P, * p_)+E4(B_ By) whole boundary 
PS,0- 1,1 S, 0+ (E _ - Eq)(E9 - E4)(E+ -E_) a,c,e+b,d,f 
A, 1+ 0,1 V,1- E_(Dy- Py) + Eo(P, - P_) +E4(P_- Pp) b,d,f only 
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FIG. 3. (A): Dalitz plot of 171 
triplets from the control region 
(820 <M3< 900); (B): folded con- 
trol region plot; (D): Dalitz plot 
for 191 triplets in the peak re- 
gion, 43 + 7% of which are due to 
w mesons; (C) folded peak region 
plot. T4,7_, and Tg are kinetic 






















energies of the 7*, 7~, and 7°, 
respectively. 


tum barrier (J>0) and constitutes mild evidence 
against a V matrix element (A meson). 

The two stronger remaining candidates have 
A vs S matrix elements. The dashed lines in Fig. 
3(D), as well as the two straight lines of the fold- 
ed distribution in Fig. 3(C), correspond to equal 
energies of two pions. The scalar matrix element 
(S) of Table I vanishes when any two pions have 
the same energy, and therefore would require de- 
population along these lines. This is not observed. 
An A matrix element has terms in p;p;, which 
vanish for collinear pions. The boundary of the 
plot represents collinearity, and seems indeed to 
be depopulated; although clearly more statistics 
and more detailed analysis such as investigation 
of polarization and alignment are needed. 

We conclude that the data fit the qualitative cri- 
teria for an axial vector matrix element (w me- 
son); there is reasonable evidence against both 
an A meson and a PS meson. 

The film used in this measurement was obtained 
in collaboration with J. Button, P. Eberhard, 

G. Kalbfleisch, J. Lannutti, G. Lynch, and N. H. 
Xuong; and this experiment would not have been 
possible without their help. It is a pleasure to 
thank Professor Murray Gell-Mann for his theo- 
retical discussions. We wish to acknowledge the 
active participation of C. Tate, L. Champomier, 
A. Hussain, C. Rinfleisch, and F. Richards in the 
final stages of this experiment. 
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TOTAL CROSS-SECTION MEASUREMENTS OF K*-p AND K*-n INTERACTIONS 
IN THE MOMENTUM REGION 0.77 TO 2.83 Bev/c* 


V. Cook, D. Keefe, L. T. Kerth, P. G. Murphy,? W. A. Wenzel, and T. F. Zipf 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received July 31, 1961) 


Total K*-p and K*-n cross-section measure- 
ments have been made previously up to 0.8 Bev/c 
with emulsions,’ counters,” and bubble chambers,° 
in the region 0.8 to 2.4 Bev/c by the M.I.T. group,* 
and at momenta above 2.7 Bev/c by the CERN ® 
and Dubna® groups. From the data of the M.IL.T., 
Dubna, and CERN experiments, it is difficult to 
arrive at a satisfactory description of the Kt-p 
cross section, o,, in the region of a few Bev/c. 
At the highest momentum measured by the M.I.T. 
group (2.4 Bev/c), 0» was found to be 12.9+1.0 
mb, whereas at 2.9 Bev/c, the CERN group found 
Op = 24.5+ 2.5 mb. The present experiment was 
undertaken partly for the purpose of resolving this 
difficulty. 

The experimental technique was similar to that 
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used by Cook et al. to measure K~~p total cross 
sections.” The beam is illustrated in Fig. 1. A 
feature of the present beam not encountered in the 
K~ beam described in reference 7 was that the 
high proton counting rate (approx 2 x10°/sec) could 
cause accidental background. This source of ac- 
cidentals was eliminated with an anticoincidence 
circuit designed to reject any K* meson that was 
accompanied by another beam particle within a 
time of +50 musec. 

The hydrogen-deuterium target was 4 ft long, 
6 in. in diameter, made of 0.007-in. stainless 
steel, and was vacuum insulated. Two transmis- 
sion counters, T, and T,, were used. The 7, coun- 
ter was a 12-in, x12-in. xj-in. scintillator, and 7, 
was a circular scintillator 9 in. in diameter and } 


\ Proton-beam 


Z . trajectory 
a ‘ 
\ 
‘\ 
Z *% Internal 
Z ‘ target 


FIG, 1. Arrangement of 
the variable-momentum K* 
beam. M;,,M,, andM 
are bending magnets; Q, 
Q, and Qs; are magnetic 
quadrupole lenses; B,, B:, 
. B;, and By, are coincidence 

hs counters; A; and A, anti- 
™ coincidence counters; and 
C, and C, are gas Cerenkov 
counters. The transmis- 
sion of the hydrogen target 
was measured by the scin- 
tillation counters 7; and 
T2. 





Qo ga 


-— pp ps ® 


a a ee ee ee 


fmth) 


oo... 





ould 
ce 


aS 


t of 
n Kt 


Me 


Bo, 
nce 
i- 
ind 
nkov 
is- 
rget 
cin- 








VoLUME 7, NUMBER 5 


PHYSICAL REVIEW LETTERS 





SEPTEMBER 1, 1961 





Table I. Total cross sections: K* - Pp» K* -d, and K*-n. The corrections that have been included are described 
in the text. The quoted errors include all known random uncertainties; in addition, there may be a systematic er- 
ror of—at most— +0.3 mb, which would translate all values together. The og, has not been corrected for Pauli 
principle effects in the deuteron. The AQ, is the solid angle subtended by 7, at the center of the hydrogen target. 





Momentum 0.77 0.97 1.17 1.30 


1.44 1.69 1.97 2.26 2.55 2.83 





(Bev/c) 
oa. (mb) 13.0 15.4 18.1 17.9 18.1 17.5 16.9 197.1 17.1 16.7 
p +0.9 +0.6 +0.6 +0.9 +0.7 +0.6 +0.4 +0.5 +0.6 +0.5 
g ,(mb) 27.3 32.4 35.4 35.6 35.4 35.1 34.6 33.9 33.4 33.4 
+0.6 +0.5 +0.5 +0.6 +0.5 +0.5 +0.5 +0.6 +0.5 +0.7 
a, (mb) 15.5 17.8 18.2 18.5 18.1 18.5 18.6 17.7 17.1 17.5 
+1,1 +0.8 +0.8 #1.1 +0.9 +0.8 +0.6 +0.8 +0.8 +0.9 
4Q,(msr) 40.0 22.5 19.8 17.8 15.9 13.6 11.9 10.0 10.0 8.6 





in. thick. The distance of T, and T, from the tar- 
get was varied between 5 ft and 10 ft, depending 
on the momentum, to minimize the corrections 
for Coulomb scattering, decay in flight, etc. Sev- 
eral alternate “target full” and “target empty” 
runs were made at each momentum. The results 
are shown in Table I and are plotted in Figs. 2(a) 
and 2(b), together with data from other experi- 
ments. 

Corrections have been included for the following: 
(a) There is a change in the decay-in-flight rate 
due to energy loss in the full target—this is larg- 
est at low momenta (approx 3 mb) but can be cal- 
culated reliably. (b) Forward scattering and for- 
ward recoil corrections to the hydrogen data were 
calculated using values for the 0-deg and 180-deg 
cross sections obtained by interpolation from an- 
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gular-distribution measurements at 1.0, 1.2, and 
2.0 Bev/c.® For the larger counter this correc- 
tion amounted, at most, to 0.4 mb. The forward 
scattering corrections to the deuterium data were 
calculated using optical theorem values for the 
forward scattering cross sections. For the larger 
counter this correction varied from 0.5 mb toa 
maximum of 1.3 mb at the highest momentum. 

(c) Multiple scattering in the target introduced 

a correction to the smaller counter results at 
770 Mev/c (0.7 mb) and 970 Mev/c (0.1 mb). (d) 
The Glauber screening correction,® used to de- 
rive o, from og -%, amounted to about 5% at all 
energies. In arriving at 0, from (0, -O%); a fur- 
ther correction should be applied for the effect of 
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FIG, 2. The K+ total cross sections as a function of momentum: (a) K*-p cross section, 0», and (b) K ton 
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the Pauli exclusion principle in suppressing for- 
ward charge-exchange scattering in deuterium. 
Until relevant data becomes available this cor- 
rection cannot be estimated accurately, and there- 
fore no correction was included. 

There is a small systematic difference between 
the cross sections measured by the counters T, 
and T,, the former giving a smaller value at each 
energy. The results given in Table I are averages 
of the corrected cross sections obtained from the 
two counters, and may therefore contain a system- 
atic error of— at most-— +0.3 mb, affecting all the 
values equally. 

Both op and o, appear to remain remarkably 
constant at about 17 to 18 mb in the momentum 
range 1 to 3 Bev/c, implying that the J=0 and 
J=1 interactions are of approximately equal 
strength. The present measurements give a val- 
ue of op above that for the highest energy point of 
the M.I,T. group, below that for the lowest energy 
point of the CERN group, and in good agreement 
with the points of the Dubna group. At lower en- 
ergies the present results are in good agreement 
with those of the M.I.T. group. The K*-p cross 
section at 0.77 Bev/c is significantly below the 
value at higher energies, in agreement with re- 
cent bubble chamber measurements.’° 

The K* -p interaction up to 0.8 Bev/c has been 
explained almost exclusively in terms of an S- 
wave interaction,’ but it is clear that the approx- 
imately constant value of op persisting to higher 
energies cannot be due to /=0 contributions alone, 
since 47x? is 18 mb at 1.0 Bev/c and 10 mb at 1.5 
Bev/c. The results of reference 8 show that there 
is a rapidly developing anisotropy in the elastic 
angular distribution at 1 Bev/c and above, indi- 
cating that higher angular momentum states con- 
tribute strongly to the elastic scattering at these 
energies. Furthermore, the K*-p interaction 
shows an increasing contribution from inelastic 
processes above 1 Bev/c. 

We wish to thank Dr. E. J. Lofgren and the Bev- 
atron Operations Group. We are grateful to Wil- 
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and to C. Burton, R. Duncan, D. Linn, and G, 
Schnurmacher who contributed to the construction 
and assembly of the experimental apparatus. 
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ANTIPROTON-PROTON AND PROTON-PROTON TOTAL CROSS SECTIONS FROM 4 TO 20 Bev/c* 


S. J. Lindenbaum, W. A. Love, J. A. Niederer, S. Ozaki, J. J. Russell, and L. C. L. Yuan 
Brookhaven National Laboratory, Upton, New York 
(Received August 16, 1961) 


Previous experiments have observed and com- 
pared the behavior of the p-p total cross section’”? 
and the p-p total cross section** from low inci- 
dent momenta (~3 Bev/c) to 10 Bev/c. Although 
the p-p cross section was slowly decreasing with 
increasing momentum from =44 mb to 40 mb, the 
p-p cross section fell rapidly from =80 mb at 3 
Bev/c to =50 mb at 10 Bev/c. The detailed be- 
havior of the p-p total cross section at even 
higher energies and comparison to the p-p total 
cross section was thought to be of considerable 
general interest, and also of specific interest 
due to a theorem of Pomeranchuk® which can be 
stated as follows: If one assumes that the total 
cross section of a particle 8 for a particle a be- 
comes constant beyond an incident energy € such 
that 


(1) 


= j > 
otal +8)=const for energies >« v 
otal +8)=const for energies >€,, (2) 


where @ is the antiparticle of a, then one can 
demonstrate on the basis of the forward disper- 
sion relations that 


otal +B) = tai @ +8) for energies >€, (3) 


where €>both €, and €,. Although the basic as- 
sumption of constant total cross sections at high 
energy which is made by Pomeranchuk’s theorem 
represents reasonable physical intuition, there is 
of course no compelling reason to believe that 
this is necessarily so, and of course no good 
reason to believe that, even if an € exists, it is 





in the 10-20 Bev energy region. 

The p-p total cross section has previously been 
shown to be constant within the errors (~5%) from 
10 to 28 Bev/c and appears to satisfy the basic 
assumption of Pomeranchuk’s theorem provided 
it remains constant at higher energies. As pre- 
viously mentioned, the p-p total cross section 
was still decreasing up to 10.5 Bev/c, the highest 
momentum previously reported, where ototal(P-p) 
- Otota] (b-p) = 10+3 mb, so that Pomeranchuk’s 
assumption was not yet satisfied. 

Using the Brookhaven 33-Bev alternating gradi- 
ent synchrotron the present investigation has de- 
termined the p-p and p-p total cross sections at 
~1-Bev/c intervals (except for two points) from 
4 Bev/c to 20.3 Bev/c. It is found that Pomer- 
anchuk’s theorem is still not verified. 

The experimental arrangement used for mo- 
menta above 9 Bev/c is shown in Fig. 1. Except 
for small corrections to the angle due to the syn- 
chrotron fringing field, the beam consists of par- 
ticles emitted from the synchrotron target at 4.5°. 
The magnets Q,, Q,, B, Q,, and Q, focus a beam 
of the desired momentum on scintillator S,, 1 in. 
in diameter. The momentum of the system is de- 
fined to ~+2.5% by the proton beam size on the 
synchrotron target and S,; the absolute scale of 
the momentum is uncertain by 2%. The particle 
compositions of these beams are given in a pre- 
vious publication.® Typical p/n~ ratios were ~1% 
at low momenta, dropping to 0.2% at 20.3 Bev/c. 
Particles in the beam are identified by the focus- 
sing gas Cerenkov counter C. A schematic dia- 
gram of the optics of the counter is shown in Fig. 





° 100 
beeeeesesd 
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FIG. 1. Experimental arrangement. T, internal target; Qi, Focussing magnets; C, collimator; B, analyzing 
magnet; S,-3, beam-defining scintillators; C, velocity-selecting Cerenkov counter; H,, liquid hydrogen target; 


S4-), scintillators for transmission measurement, 
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2(a). In Fig. 2(b) is shown a curve of the ratio 
S,S,S,CA /S,S,S, in the region of pressure corre- 
sponding to antiprotons at 16 Bev/c. As can be 
seen, there is a small background (~4%) at pres- 
sures at which the counter should not respond. It 
has been verified that these counts are due to in- 
complete rejection of pions and K mesons, and a 
correction of ~2% was made to the cross section 
for this effect. To measure the cross section al- 
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FIG, 2. (a) Velocity-selection Cerenkov counter. 
The solid ray is a typical ray from a particle with the 
desired velocity, and the broken ray is that from a par- 
ticle with higher velocity. The first is counted in chan- 
nel C, the second in channel A, (b) Cerenkov counter 
response in the region where 16-Bev/c antiprotons are 
expected to be detected. 


186 


ternate runs are made with a 10-foot long liquid- 
H,-filled target and two Al plates which simulated 
the empty target. By requiring, in addition to the 
(S,S,S,CA) identification, a coincidence with one 
of the counters S,-S,, ranging from 4 to 12 in. 

in diameter (see Fig. 1), simultaneous measure- 
ments are made of the transmission of the hydro- 
gen for six different geometries. A seventh chan- 
nel measured accidental coincidences with the 12- 
in. counter delayed. The transmissions thus meas- 
ured, of course, include particles scattered - 
through small angles as well as secondaries pro- 
duced in inelastic collisions in the target. 

By plotting the transmissions of the various 
counters S, -S, against their respective areas 
and extrapolating to zero area, one obtains the 
proper transmission which should result from an 
experiment having an ideal “good geometry.” 

Cross sections at 4-10 Bev/c were obtained 
with an arrangement similar to that shown in Fig. 
1 except that the beam used leaves the synchro- 
tron target at 9° and the angle of bending in mag- 
net B is 5.5°. The momentum resolution of this 
beam is +1.5%. The Cerenkov counter described 
above was used to measure 6-10 Bev/c while at 
4, 5, and 6 Bev/c a pair of smaller counters ca- 
pable of withstanding higher gas pressures were 
used in tandem in order to obtain a high rejection 
ratio for pions. 

At the lower momenta the effects of multiple 
Coulomb scattering in the hydrogen target are 
not negligible. Sternheimer’s method’ of cor- 
recting for this effect was used. The smallest 
counter at low momenta requires an appreciable 
correction and was used only as a monitor on the 
beam size during runs with the hydrogen removed. 

The cross sections shown in Fig. 3 were ob- 
tained by plotting the ratio of transmission with 
hydrogen in to transmission with hydrogen out 
versus counter area; each transmission was cor- 
rected for accidental coincidences first. These 
were small for the p-p but were as much as a 
few percent for some cases of p-p at low momen- 
ta where the proton beam rates were high. The 
corrected ratios, as is expected, lie on a linear 
plot since the angles subtended by the counters 
S,-S, are small enough to be within the isotropic 
portion of the diffraction scattering peak. The 
smallest counter used in the extrapolation con- 
tains all but a negligible number of the singly 
Coulomb-scattered particles and hence no cor- 
rection is required for this effect. 

In Fig. 3 the errors shown are a combination 
of statistical errors and an estimated error on 
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FIG. 3. Total p-p and p-p cross sections. 


the extrapolation to zero solid angle, an error 
associated with the multiple Coulomb scattering 
correction, and background errors. 

It is clear that our results generally agree with 
the early measurements? of total P-p) up to 10.5 
Bev/c, and with the earlier ototai(>-p) up to 28 
Bev/c* which used a C-CH, difference instead of 
a liquid hydrogen target. There is no evidence 
from our measurements for any appreciable 
structure in either ototal(p-p) or ototal(p-p). In 
fact, we also find that ototai(p-p) is constant 
within the errors from 10-12 Bev/c to 20 Bev/c 
with a value o¢otai(>-p) = 39.541 mb. On the 
other hand, ototal(p-p) appears to decrease 
monotonically from 10 Bev/c [where ototal(p-p) 
=5844 mb] to 20.3 Bev/c [where ototai(P-p) = 48 
+4 mb]. 

Amati, Fierz, and Glaser have shown® that in 
order to assure constancy of the cross sections 
at high energies and convergence of the integral 


t —— co 
sf on(8) OB ae - 





in which the particle cross section is o*t(E) and 
the antiparticle cross section is o (E), one must 
require that 


[o7(E) -o*(E)]-0 (5) 


faster than 1/ln E as E~~, The convergence of 
the above integral is necessary to give meaning- 
ful results for many dispersion relation treat- 
ments. 

If we assume that the p-p and p-p cross sections 
above 10 Bev/c have reached their limiting ener- 
gy-dependence region, and if we assume the de- 
pendence to be [o~(£) - o*(E)]~«1/ln E, then a fit 
of our data implies that o~(£) - o*t(E) may still be 
greater than 1 mb even for energies of several 
hundred Bev. Even if o~(E) -o*(E) were to ap- 
proach zero as fast as linearly with energy, it 
would appear from our data that the difference 
would not become zero before 35 Bev. 

Von Dardel et al.” showed that o(K~-p) was sub- 
stantially larger than o(K*-p) up to 8 Bev/c, where 
Stotal (K ~-P) - ototal(K*-p) = 542 mb. The results 
of their experiments as to whether o¢ota)("~-p) 
and ototal(™*-p) were the same or different up to 
10 Bev/c were inconclusive due to possible exper - 
imental uncertainties. It was also not clear in 
these experiments that the cross sections at the 
higher energies were energy independent, but at 
least their energy variation appeared to be slow. 
Therefore, to date there has not been a clear dem- 
onstration of even a single case where Pomeran- 
chuk’s theorem was apparently applicable and 
verified. 

The authors wish to acknowledge the coopera- 
tion of the AGS staff and the operating crew dur- 
ing these experiments. The authors also wish to 
acknowledge the assistance of J. Elias and C. Sin- 
clair during these experiments. 
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K’ -p INTERACTION AT 455 Mev* 
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Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 


and 


William Slater, Donald M. Stork, and Harold K. Ticho 
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(Received July 31, 1961) 


We have undertaken a systematic study of the 
interaction of positive K mesons with hydrogen 
and deuterium in the energy interval from 0 to 
455 Mev. Some preliminary results have already 
been reported.’ In this note we present the re- 
sults obtained for the elastic and inelastic Kt-p 
interaction at 455 Mev. 

Previous investigations of K* interactions in 
nuclear emulsion? and propane® have yielded 
measurements of differential and total cross 
sections for the process K*+p+K*+p in the 
energy interval 40 to 300 Mev. The differential 
cross section at 225 Mev,* as well as the total 
cross sections in the range 175 Mev< Ex <8 Bev, 
have been measured by counter techniques.*”® 
The features of the K*-p scattering from 80 to 
300 Mev are that (a) the total cross section is 
approximately 14 mb, varying little, if at all, 
with energy; and (b) the angular distribution is 
isotropic. 

The Lawrence Radiation Laboratory 15-inch 
hydrogen bubble chamber was exposed to a sep- 
arated beam of K+ mesons produced by the 6- 
Bev circulating protons of the Bevatron (Fig. 1). 
The system was designed for a momentum of 
645 Mev/c. With adjustment of the magnet pa- 
rameters it was possible to obtain the higher 
momentum of 810 Mev/c (Tx =455 Mev). A 
mass-resolution curve at 810 Mev/c for the 
separation system is shown in Fig. 2. The back- 
ground of light particles (pions, muons, and elec- 
trons) was approximately 10%; the pion compo- 
nent is analyzed in more detail below. 

The initial sample of events, chosen to satisfy 
geometric and incident-momentum criteria, con- 
tained both inelastic and elastic K* interactions 
and also a background of 7* interactions. Kine- 


188 











Scale 
Oo 5 1 Feet 
-—+— 


@ 5 
“—————-15-inch bubble chamber 


FIG. 1. Layout of the separated K* beam, The beam 
design was similar to a separated K~ beam designed 
earlier [P, Eberhard, M. Good, and H. K. Ticho, 
Rev. Sci. Instr. 31, 1054 (1960)], and is described in 
detail in G. Goldhaber, S. Goldhaber, J. Kadyk, T. 
Stubbs, D. H. Stork, and H. K, Ticho, Lawrence 
Radiation Laboratory Internal Report Bev-483, Febru- 
ary 26, 1960 (unpublished). The K* beam from the 
target (7) is focused by the quadrupole Q, onto slit S;. 
The momentum selection is effected by bending magnet 
BM,, and the subsequent mass separation by the 
crossed electric and magnetic field in spectrometer 
Sp,. The second state is essentially a mirror image of 
the first. The steering magnet SM was introduced for 
additional freedom in the horizontal plane. Cpoyjz and 
Cyvert re horizontal and vertical collimators, respec- 
tively. 
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FIG, 2. Mass analysis of particles emerging from 
slot S; in Fig. 1. This curve was obtained by setting 
spectrometer Sp, to transmit K mesons and varying the 
magnetic field in spectrometer Sp,. One thus obtains 
amass analysis of particles leaving slit S,;. The final 
operating conditions for Sp, are indicated by the arrow. 


matical fitting procedure (with further evidence 
from estimates of bubble density of the secon- 
dary tracks) provides means of separating the 
elastic from inelastic interactions with essen- 
tially complete certainty. The pion elastic scat- 
terings are distinguishable from K* elastic scat- 
terings for cos6,°-™-<0.4. To determine the 
total number of pion scatterings, we use the 7*+-p 
angular distribution’ to evaluate the number of 
such events in the remaining angular interval 0.4 
<cos6,°"™-< 0.96. This number (i.e., 31 events) 
is then subtracted from the group consistent with 
K*-p elastic scatterings. The observed number 
of 7* inelastic interactions (i.e., 29 events) is in 
agreement with the sum of observed and inferred 
numbers of 1* elastic interactions. The contami- 
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nation of pion inelastic scatterings in the sample 
of K*-p elastic scatterings is thus negligible. We 
accepted only those scatterings with cos6,°-™- 

< 0.96 in order to preclude any effects of low ef- 
ficiency for detecting small-angle scattering. We 
find, then, a total of 1330 elastic K*-p scatterings 
satisfying the above conditions. The total K* path 
length was determined by three independent meth- 
ods: 

(a) From K* decays into three charged secon- 
daries (316 decays) and the known branching ra- 
tio,® b, = 0.0614 0.002 for these decay modes; this 
yields a total path length of (3.124 0.25) x10° cm. 

(b) From Kt decays into one charged secondary 
with projected angle 6),,>27.5 deg. This cutoff 
angle was introduced in order to avoid possible 
confusion with K*-p scatterings without observa- 
ble recoil. The decays included in this sample 
correspond to a fraction 6, =0.29+ 0.01 of all K* 
decays. A total path length of (2.93 + 0.14) x10® 
cm is obtained by this method. 

(c) From a direct count of tracks passing 
through the chamber. After corrections for light- 
particle contamination, decays, and interactions, 
we obtain a total path length of® (2.96 + 0.16) x10° 
cm. The weighted average of the path lengths ob- 
tained by the above three methods is Ltotg) = (2.97 
+ 0.09) x10® cm. Correcting for the scanning ef- 
ficiency for elastic scattering, €,=0.997, and 
decays, €g=0.994, and extrapolating the angular 
distribution (considered as flat) to cos@,. m, =1.0, 
we obtain the total elastic-scattering cross sec- 
tion, at 455+ 5 Mev, og; =13.0+0.7 mb. This 
cross section contains some Coulomb effects be- 
low cos6x°*™: =0.96, and is to be compared with 
the purely nuclear cross sections as deduced from 
the phase-shift analyses below. The error in- 
cludes the uncertainties in the path-length deter- 
minations and statistical error in the number of 
scatterings. For the inelastic K*-p interactions, 
discussed below, we obtain a cross section On) 
=1.0+0.2 mb. 

The differential cross section plotted in Fig. 3 
shows only a small angular dependence. We ana- 
lyzed the data in terms of s- and p-wave scatter- 
ing for the elastic and inelastic interaction. 

We can thus write the differential cross section 
as 


l+ia 246. 


2 2 
“1+; Ty (Ma? 1 + 2m5e Sas _ 3) cosé 
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FIG. 3. The K*+H elastic dif- 
ferential cross section at 455 
Mev. The results correspond 
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to 1320 scattering events. The 
curves are computed from the 
various “best fit” phase shifts 
as given in Table I. Sets B’, 
Cc”, and B*C* give essentially 
identical differential cross-sec- 
tion curves. For clarity, only 
— Set A™ and Set B’ and C~ are 
shown in the main figure. The 
inset shows the small-angle be- 
— havior of all the phase-shift 
solutions. 
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and the inelastic cross section as 


Oinel = mA = n,”) + (1 =") + 2(1 -™s°) ]. (2) 
Here a =e? Mv re) 6 is the c.m. scattering an- 
gle, ik the c.m. momentum, and vre] the relative 

velocity’; 6,, 5,,, and 5,, are the Sy», Py2, and 
by2 T=1 phase shifts, respectively, and 7,, 1, 
and 7,, correspond to the imaginary part of these 
phase shifts. For simplicity we have assumed 
that the inelastic scattering occurs principally 

in one of the three scattering amplitudes. Solu- 
tions for the phase shifts were obtained by setting 
two of the absorbitive amplitudes 7, equal to unity 
and obtaining the third one from (2). The solutions 
are rather insensitive as to which phase shift was 
chosen as complex. In Table I and Fig. 3 we give 
the solutions corresponding to n, =0.92+0.2 and 
Tha =™s =1. 

There are three sets of phase-shift solutions. 
Set A: A dominant s-wave solution. The sign 
of the s,,., phase shift can be seen to be most prob- 
ably negative in agreement with earlier results 

at lower energy.**™ 

Set B: A dominant p,, solution which is the 
Minami ambiguity corresponding to Set A. A 
unique determination of the sign of 5,, from our 
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data is not possible because Coulomb interfer- 
ence here occurs at smaller angles than for solu- 
tion A. 

Set C: A combination of p,,. and py, amplitudes 
such as to reproduce near isotropy with an am- 
biguity in sign. This is the Fermi-Yang ambigu- 
ity corresponding to solution B. If we consider 
the evidence for a repulsive (i.e., positive) nu- 
clear potential from the emulsion data?*™ whose 
largest contribution comes from the forward 
scattering amplitude in the T =1 state; viz., V 
~ -Re[0.75f,(0) + 0.257,(0)],we can infer the sign 
of the dominant phase shifts. For reasonable 
values of f,(0) this would rule out solutions A’, 
B*, and C’, 

Isotropy and little variation of the scattering 
cross section dominate the K*-p interaction 
throughout the energy interval up to 455 Mev. 

To ascribe the scattering to predominant re- 
pulsive s-wave interaction even at 455 Mev 

(Set A~) does imply an anomalously low p-wave 
interaction. A p,,. solution (Minami ambiguity) 
can clearly fit an isotropic distribution at any 
energy. The near constancy of the cross section, 
however, over the large energy interval makes 
a dominant p,, (Set B) or py - py, mixture (Set C) 
solution rather unlikely. At this point we would 
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Table I. Phase shifts for K’ -nucleon scattering at 455 Mev in the T=1 state.” 








Phase shifts 





(deg) Probability “a 

Solution 5, 511 5:3 from a x’ fit (mb) 
A~ -47 +1 0.5+4.5 1.5+2.5 0.15 12.2+0.4 
AY 49 +1 - 0.544 0+2 0.01 13.0+0.4 
B™ 4.542 -45.5 +1 - 2.541 0.92 12.6 +0.5 
Bt - 1,541.5 46.5 +1 4 #1 0.73 13.0 +0.5 
c 4.542 14.5 +1.5 -28 +1 0.92 12.5 +0.5 
c* - 1.542 -13.5 +1.5 29 +1 0.73 13.0+0.5 





“The solutions given here are computed for 7;=0.92, m4;=13=1. If we take 744=0.92 and n;=43=11, for ex- 


ample, Solution A’ becomes 6; =-45°, 64,=3.5°, 643=1.5°. 


This is the nuclear elastic cross section computed from the respective phase shifts leaving out the Coulomb 
terms. The errors reflect the errors on the respective phase shifts. 


like to emphasize that we have not explored com- 
binations of s, p, and d waves, which will, of 
course, also reproduce the K*-p scattering proc- 
ess. As is well known from the proton-proton in- 
teraction, certain combinations of several angu- 
lar momentum states can reproduce isotropic 

and energy-independent differential cross sec- 
tions over appreciable energy intervals. 

As can be seen from the inset in Fig. 3, a pre- 
cise measurement of the scattering at small an- 
gles can distinguish between dominant s and dom- 
inant p solutions because of the difference in the 
respective Coulomb interference. Similarly, po- 
larization measurements of the recoil proton 
could determine the presence of a mixture of 
by2 and py, amplitudes or possibly higher partial 
waves. Scattering in the s,,. and in the pure p,,. 
states, respectively, does not give rise to polar- 
ization. 

Inelastic interactions of positive K mesons with 
single pion production can proceed via three pos- 
sible channels. Among the 102 inelastic scatter - 
ings recorded in the chamber, we observed ex- 
amples of all three modes of pion production. 
Table II summarizes the results. 

Those examples of Reaction I with subsequent 
K,° decay are readily identifiable; 25 were ob- 
served. For a K,° branching of § into charged 
pions, these events should represent 3 of all K° 
mesons produced in Reaction I. If all ambiguous 
inelastic scatterings belong to channel I(b), the 
observed ratio is consistent with that expected. 

In any event it is clear that Reaction I dominates 
strongly. It is interesting to note that if the K 


and 7 mesons were produced in a total isotopic 
spin state of T= 4, the ratio of the rates of Reac- 
tions I, I, and II would be 2:1:0, while produc- 
tion of the 7 meson and nucleon in the T = 3 state 
would yield at ratio 9:2:1. The data are sugges- 
tive of a dynamical effect which may be due to 

an enhancement in the production in one of these 
isotopic spin states. 

We wish to thank Professor Luis W. Alvarez 
and many members of his group for making the 
15-inch bubble chamber and analyzing facilities 
available to us. We are very grateful for the 
tireless efforts of the bubble chamber crew and 
the Bevatron crew as well as our own scanning 


Table I. Pion production in K*- p collisions. 








Number 
Channel of events 
I (a) K*+p—K,°+n* +p 25 
\ + @ 
T +7 
(b) a’ +p— Ken +p 35 
K,° or neutral decay of K,° 
IK +p—-K*+n+p 24 
Wl K’+p—-K*+n*+n 8 
Ambiguous—I(b) or II " 10 
Total 102 











*Included in this number are nine events also con- 
sistent with 1-p inelastic scatterings. 
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and measuring group, without whose assistance 
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1-1 RESONANCE IN 1~-p INTERACTIONS AT 1.25 Bev* 


E. Pickup,f D. K. Robinson, and E. O. Salant 


Brookhaven National Laboratory, Upton, New York 
(Received August 14, 1961) 


It has become apparent that 7-7 interactions 
are important in 7-p collisions in the Bev region. 
Several experiments on single-pion production’ 
have shown that there is an excess of nucleons at 
low laboratory kinetic energies over that which 
would be expected on a statistical or isobar theo- 
ry. Goebel’ has pointed out that such an excess 
would be indicative of a 7-7 interaction. Fraser 
and Fulco® showed that a 7-7 resonance in an 
I=1, J=1 state could explain features of elec- 
tron-nucleon scattering, and suggested a reso- 
nance at w*~ 10u?, where w is the total energy of 
the two pions in the 7-7 rest frame, and yp is the 
pion rest mass (w=Q+2u). Later calculations® 
suggested that the resonance was at w*~ 22y?. 
Indications of a 7-7 resonance in 7-p interac- 
tions have been obtained in several experi- 
ments,°’*#”” but the results were limited by low 
incident pion energies, or by statistics. With the 


192 


higher statistics available in the present work, 
strong evidence is obtained for a 7-7 resonance 
in an J=1 state. 

As part of a study of 7-p interactions, using 
the bubble chamber technique, we have meas- 
ured 4000 7~-p events at 1.25-Bev incident pion 
energy. 968 events of the type 


17 +pon-+mten, (1) 
and 566 of the type 
17 +p2m-+7°+p, (2) 


have been identified using the GUTS kinematic 
fitting program, and ionization density measure- 
ments. 

Comparison of the pion momentum spectra with 
the predictions of the extended isobar model** 
indicate that pion production through isobar for- 
mation is not the dominant process at this energy. 
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The nucleons from both reactions are produced 
predominantly backward in the c.m. system, and 


thus are peaked at low laboratory kinetic energies. 


Figures 1(a) and (b) show the nucleon laboratory 
kinetic energy, 7, for reactions (1) and (2), re- 
spectively. The curves in Fig. 1 were computed 
assuming a constant 7-7 cross section, 07-7, 
from the equation given by Chew and Low”"*: 


aA" Ow" 


f? A?/u? w 
A? + - py? 2a (47+ 2”) WL 


where A is the four-momentum transfer to the 
nucleon (A?=2MT), f? is the renormalized pion 
coupling constant and q,; is the laboratory mo- 
mentum of the incident pion. The theoretical 
curves are normalized to the total numbers of 
events in the histograms. The theoretical and 
experimental distributions seem to agree at low 
momentum transfers, but the histograms are 
more sharply peaked than would be expected.*® 


(4w?-u?)o_, (3) 
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FIG. 1. Histogram of nucleon laboratory kinetic en- 
ergies. The curve was calculated from the Chew-Low 
equation, for constant 0,_,. 


There is also a group of high-momentum-trans- 
fer events, not expected on a one-pion exchange 
model. 

Figure 2 shows Q values for 7~-1* and 1~-1° 
systems, and also the prediction of the statis- 
tical model. The isobar prediction is similar 
to the latter. Both histograms exhibit well- 
defined peaks at the same @ value. The agree- 
ment in position of the two @ peaks suggests that 
the same 7-7 resonance state is dominant in both 
reactions. It follows, therefore, that isotopic 
spin J=0 is excluded for this state. 

The branching ratio for reactions (1) and (2), 
R=0(2)/o (1) =0.585+ 0.031, and the ratio, Ry_,, 
for events in the Q,_, peaks (400 Mev <Q,_, 
<550 Mev), is 0.504 0.04. The one-pion exchange 
model predicts that Rz_, = 4 for a 7-7 interaction 
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FIG. 2. Q-value distributions for 7-1 system for all 
events, for (a) m-+p— 1-+nt++n, (b) m-+p—1-+7°+p; 
the dashed curves are from statistical theory. 


193 











































VOLUME 7, NUMBER 5 


PHYSICAL REVIEW LETTERS 


SEPTEMBER 1, 1961 





in a pure /=1 state, and § for J=2. The model 

is more likely to be valid if we select events with 
small A”. Restricting the selection of events in 
the Q peak to those with A?<10u?, or T<104 Mev, 
we obtain R,_, =0.42+0.06. The data are in 
agreement with the predictions of the model for 
I=1. This state is further suggested by the Yale 
m+-p experiment at 1.26 Bev, ** by comparing the 
n*-n° and m+-2* Q distributions. 

The 1-7 cross section, as a function of w?, was 
calculated by integrating Eq. (3) up to A? =10y’, 
assuming this equation to be valid in the physical 
region. Comparison of the experimental and the- 
oretical distributions in Fig. 1 indicates that this 
may be a reliable procedure for small A?.*® Fig- 
ure 3 shows 07-7 obtained in this way, for the 
combined events from reactions (1) and (2). The 
background on the low-energy side of the reso- 
nance is associated with both reactions. This 
would seem to eliminate a pure /=0 low-energy 
effect. We obtain o7., =95 mb at the peak, w? 
=29u? (Q=475+10 Mev). The full width at half 
maximum is 130 Mev. These values, and the 
values determined at 1.74 Bev by Erwin et al.” 
agree within statistics. This agreement between 
results at two quite different energies implies 
that there is, indeed, a resonance, 7*, in the 


1-7 system. 

The decay angle, 5, between the direction of 
the 7~ in the rest frame of 7* and the c.m. direc- 
tion of ™*, was calculated. In Figs. 4(a) and (b) 
the distributions of events in the Q peaks show 
backward-forward symmetry, but a marked de- 
parture from isotropy. Restricting the selection 
to events with small A? shows that the anisotropy 
is mainly associated with such events. For events 
outside the Q peak, the asymmetries observed 
in both reactions (1) and (2) are consistent with 
some (§, 3) 7-p isobar formation. 

To investigate angular momentum states in- 
volved in the 7* resonance, the angle, a, between 
the direction of the 7~ in the rest frame of 7* and 
the direction of the incident pion, was calculated. 
Events were selected in which 7* was emitted 
closely parallel or antiparallel to the incident pion 
direction ([cos@|>0.9).47 Provided that the Adair 
analysis” holds for this selection, and for a short- 
lived resonant state, the 7~ angular distribution 
is expected to be a relatively simple function of 
the angular momentum state J (a+b cosa for 
J=1). The distributions are shown in Figs. 4(c) 
and (d). That in 4(c), at least, is suggestive of 
J=1, rather than larger J values which would give 
higher powers of cos*a. The distribution in 4(c) 
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FIG. 4. Angular distributions for the 7-7 system. 

(a) and (b) show decay angles 6 of 1- in 7* rest frame, 
relative to c.m. direction of 7*: (1) All events (dotted); 
(2) events with 400 Mev< Q,_, = 550 Mev (dashed); (3) 
events from (2) with A?< 10? (solid). (c) and (d) show 
decay angles a of m~ in 7* rest frame, relative to inci- 
dent pion direction, for events with |cos@|>0.9, where 
6=c.m. angle of 7*. 


shows some asymmetry (backward/forward ratio 
=43/76), which may indicate interference from 
another state. 

We wish to thank the BNL bubble chamber 
group, the Yale group who set up the beam, and 
our efficient team of scanners. 
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ONE-PION EXCHANGE IN p-p COLLISIONS AT 2 Bev* 


W. J. FickingerT 
Yale University, New Haven, Connecticut and Brookhaven National Laboratory, Upton, New York 


and 


E. Pickup, D. K. Robinson, and E. O. Salant 
Brookhaven National Laboratory, Upton, New York 
(Received August 14, 1961) 


Recently several authors'~* have proposed the 
use of the one-pion exchange model to describe 
the mechanism of pion production in p-p colli- 
sions. It is assumed that the dominant process is 
the exchange of a single pion between nucleons, 
with scattering of the virtual pion on one nucleon. 
The model then allows the calculation of the nu- 
cleon angular and momentum distributions, and 
predicts the absolute pion production cross sec- 
tions, using the known 7-) elastic scattering cross 
sections and the renormalized pion coupling con- 
stant. Experimental results of Batson et al.° and 
Fowler et al.* have been compared with the predic- 
tions of the model by Kobayashi, Selleri,? and Ii- 
zuka and Klein.‘ However, the limited statistics 
in these experiments do not provide an adequate 
test of the model. 

In a recent experiment at the Cosmotron, pro- 
ton-proton interactions at 2 Bev were observed 
in the 20-inch hydrogen bubble chamber. The re- 
sults of a complete analysis of this experiment 
are being prepared for publication.” From the 
measurement of about 4000 two-prong events, 
1300 events of the type 


p+p>~p+ne+nt 


were identified, using digitized equipment and the 
GUTS kinematic fitting program, together with 
ionization density measurements. The results 
have been compared with the theory presented by 
Selleri. In general the agreement between theory 
and experiment is remarkable, but there are sig- 
nificant deviations. 

The histogram in Fig. 1 shows the experimental 
c.m. momentum distribution of the neutrons. The 
dashed curve, normalized to the number of events, 
gives the prediction of the statistical theory, while 
the solid curve was computed by numerical inte- 
gration of Selleri’s equation,® using the experimen- 
tal 1*-p elastic scattering cross sections. The 
contribution from 7°-p charge exchange scatter - 
ing has been neglected in the calculation. Fora 
pure T=3/2 state, this is 1/9 the contribution 
from 7-p scattering. It should be noted that, 
since absolute cross sections are given by the 
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FIG, 1. Variation of partial cross section with neu- 
tron c.m. momentum for the reaction, The histogram 
is the experimental distribution; the dashed curve is 
from statistical theory (normalized); the solid curve is 
calculated from the equation of Selleri. 


theory, this latter curve is not normalized to the 
data. The sharp peak at 800 Mev/c in the theo- 
retical neutron momentum spectrum is a conse- 
quence of the (3/2,3/2) pion-proton resonance. 
A similar peak is also predicted by the isobar 
model of Lindenbaum and Sternheimer,’ although 
this model does not give absolute cross sections 
or production angular distributions. 

The c.m. angular distribution of the neutrons 
(folded about 90°) is shown by the histograms in 
Fig. 2. The dashed curves were computed from 
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FIG. 2. Variation of cross section with c.m. angle 
of the neutron compared with the prediction from Sel- 
leri’s equation. (a) shows the distribution for |cos6| 
20.8; (b) shows the distribution for all angles, folded 
about 90°, 


Selleri’s equation, thus giving the absolute dif- 
ferential cross sections. The model gives a rea- 
sonable fit to the data at small angles, although 
the predicted differential cross section drops off 
less rapidly than the observed distribution. At 
large angles the theoretical cross sections are 
about twice those observed experimentally. Thus, 
the total cross section predicted was 23.4 mb, 
compared with our experimentally determined 
value of 16.06+ 0.44 mb for this reaction. How- 
ever, events with large recoil angles involve 
large momentum transfers, i.e., collisions with 
small impact parameters. Consequently, the con- 
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FIG. 3. Scattering angle of the virtual pion in the 
n*-pc.m, system, (a) Events outside (3/2, 3/2) 
resonant peak; (b) events in resonant peak. 


cept of real scattering of an intermediate virtual 
pion is unlikely to be valid in such collisions. 

An independent test of this concept might be giv- 
en by considering the angular distribution of scat- 
tering of the virtual pion. If these can be consid- 
ered to behave as real pions the characteristic 
(1+3cos*q@) distribution, associated with 7-p 
scattering in the angular momentum state J =3/2, 
should be observed for those events for which the 
total kinetic energy of the pion and proton in their 
c.m. system, Qzp, is in the region of the (3/2, 3/2) 
resonance. The direction and momentum of the 
incident virtual pion were assumed to be given by 


197 








VOLUME 7, NUMBER 5 


PHYSICAL REVIEW LETTERS 


SEPTEMBER 1, 1961 





the resultant of the laboratory momenta of the 
final-state pion and proton. The angle, a, of the 
final-state pion in the 7-p center-of-mass system, 
relative to this virtual pion direction, was calcu- 
lated for each event. Of two possible rest frames 
(incident proton and target proton) for each event, 
that one was chosen in which the neutron retained 
the higher momentum, so that the virtual pion 
production was more probably associated with a 
peripheral collision. Events were selected from 
the resonant region with 120 Mev <Qzp< 200 Mev. 
An examination of the experimental 1*-p scat- 
tering data’® showed that the mean angular dis- 
tribution through this region is approximately 

(1+ 2.5 cos*q). 

The results are shown in Fig. 3. Events in the 
selected Q interval were approximately symmet- 
ric about 90° (forward/backward = 254/227). The 
folded distribution shows partial agreement with 
a (1+2.5cos*q@) distribution. Events with Q<120 
Mev give isotropy within statistics, and for those 
with Q 2 200 Mev the distribution is peaked strong- 
ly forward in qualitative agreement with 1*-p scat- 
tering above the resonance. 

The shape of the neutron momentum distribu- 
tion predicted by the one-pion exchange model is 
in good agreement with experiment. The theoreti- 
cal differential cross section is too high at large 
neutron angles, so that the total cross section is 
about 50% higher than that observed. At small 
angles (small momentum transfers) the mean dif- 
ferential cross section observed is in good agree- 
ment with theory, but the cross section changes 
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more rapidly than the theory predicts. 

We wish to thank the BNL bubble chamber group, 
our capable scanners, and Dr. Philip Duke for 
discussion and for information on 7*-> scatter- 
ing cross sections. 
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QUASI-ELASTIC PEAK IN HIGH-ENERGY NUCLEON-NUCLEON SCATTERING* 


S. D. Drell and K. HiidaT 
Institute of Theoretical Physics, Department of Physics, Stanford University, Stanford, California 
(Received August 4, 1961) 


A bump in the energy spectrum of inelastically 
scattered protons emerging from proton-nucleus 
collisions for incident energies in the range 9-25 
Bev, and for 20-60 milliradians (mr) scattering 
angles has been reported by Cocconi et al.’ at 
CERN. Characteristic feature of this bump are: 

1. The energy difference between the inelastic 
peak and the peak corresponding to elastic scat- 
tering at the same angle is =0.8 -1.1 Bev and is 
roughly independent of the scattering energy and 
angle. 

2. The height of the inelastic peak falls off 
with increasing momentum transfer as does the 
elastic cross section, in a manner suggestive 
of diffraction scattering. 

It is the purpose of this Letter to suggest a 
mechanism giving rise to a bump with these 
characteristics. This mechanism is the diffrac- 
tion scattering of a pion in the cloud of the target 
nucleon from the incident projectile nucleon. It 
leads to a simple physical picture with an appar- 
ently general application. We wish to discuss 
this here, and to give the main results of calcu- 


lations for the conditions of observation at CERN. 


A detailed report of the calculations and espe- 

cially of the relevant, and somewhat involved, 

relativistic kinematics is in preparation.? 
Consider the Feynman diagram of Fig. 1 which 


FIG, 1. Feynman dia- 


gram for scattering of a P; (A) k 
pion in the cloud of the 

target nucleon from the ~~@, Pr 
incident projectile nucle- 4 


on. The circle o on the 

nucleon line denotes the 

high-energy incident and q: (B) q 
observed final nucleon. ’ 


Pj 


2 
m 
» 


(a) Pe (b) q 








FIG, 2. Feynman diagrams for contributions neglected in the present calculation. 


introduces the notation q;, pj, and p¢ for the four- 
momenta of a target, incident, and observed final 
nucleon, respectively; 9;*=p;7= bf =m’*; fi=c=1). 
q is the recoil nucleon four-momentum and A =q; 
-q is the transferred momentum; A?<0. The 
meson produced in the collision carries off k, 
with k? =. In the following, the binding of the 
target nucleon in the target nucleus is neglected 
and the mass ratio n/m —0. 

We are interested only in a very restricted re- 
gion of scattering through a laboratory angle 6, 
<60 mr and with small energy loss € =5+t/2m 
=1 Bev«E;, where t=(; -p,)*<0 is the invari- 
ant momentum transfer between the incident and 
observed nucleons of energy E; and Ef, respec- 
tively, in the lab system, and 6=E; -Ef; €=0 
for elastic scattering. For such events the mag- 
nitude of the momentum transferred between the 
two vertices A and B can be small [lAl min’ 
~(2.4u)? for E;=25 Bev, 6,=40 mr, and €=1 
Bev] and we approximate the scattering ampli- 
tude to the pole term shown in Fig. 1. 

This approximation has been widely discussed 
recently’ and is applied here with optimism in its 
at least qualitative validity since there is an en- 
hancement factor at vertex A which we discuss 
in the following paragraph. We neglect contribu- 
tions from graphs of the types in Fig. 2. Figure 
2(a) is an exchange contribution and is unimpor- 
tant when one of the two nucleons has almost all 
of the energy. Figure 2(b) was discussed pre- 
viously* and is discarded here for two reasons: 
The momentum transfer between vertices is =1 
Bev and far from the peripheral condition of an 
almost real pion being exchanged. If, neverthe- 
less, the pole approximation is applied, its con- 
tribution to the inelastic peak is smaller than the 
observed cross section as well as smaller than 


} (n) 


(c) (d) 
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the contribution calculated from Fig. 1. (For Ej 
=16 Bev, 6,=40 mr, and €=1 Bev, it is smaller 
by a factor =4; see also footnote 10.) Feld and 
Iso® have shown that the inelastic bump occurs 

at an energy corresponding to a mean of the sec- 
ond and third resonances in pion-nucleon scatter- 
ing and suggest this isobar formation as the cause 
of the observed inelastic bump. Their argument 
is kinematical, without reference to specific 
(poor) field-theoretic approximations. However, 
it is difficult to reconcile with the absence of an 
inelastic bump at € = 300 Mev corresponding to 
the large 3-3 resonance in pion-nucleon scatter- 
ing. Figure 2(c)* can be shown’ to contribute in- 
significantly for the small values of € =1 Bev of 
interest here, primarily due to phase-space lim- 
itations. Figure 2(d) for more than a single-pion 
exchange between A and B is not calculable and 
its neglect can be defended only by measuring 
the ability of Fig. 1 in reproducing the observed 
cross sections. 

An inelastic bump at the observed € and of the 
observed magnitude emerges from Fig. 1 in the 
following way. The interaction at vertex A is 
treated as physical pion-nucleon scattering for 
small values of |A?|<M?. The spin average of 
the square of the invariant matrix element for 
the pion-nucleon scattering, 22/spins!MA|*=Xa, 
can thus be related to experimental cross sec- 
tions. 

In the rest system of the incident nucleon we 
can write 


t d t 
x -42n*(- ) "A, 1ab®1" (1) 
A S, -m* +t dQ : 





in terms of the invariant total energy s, = (p¢+k)* 
>0 and momentum transfer ¢ = (p; - p7)?< 0. We 
are interested in the experimental conditions at 
CERN which limit momentum transfers to |t| <2 
(Bev)?. On the other hand, the scattering energy, 
s,, can be much larger for inelasticities €=1 
Bev. It is this possibility that leads to the en- 
hancement at the vertex A: For very large en- 
ergies s,, the pion-nucleon differential cross 
section develops a large forward diffraction peak. 
For forward scattering, ¢=0, and we have in the 
approximation that there is mainly diffraction 
scattering at high energies: 





do, tab 1?” ~ “tot *(s,) (s, -m*)? (2) 
an aan)" om 


Here o¢o¢(S1) is the total cross section for the 
pion-nucleon cross section and is nearly con- 


200 


stant” [o¢o¢ = 25-30 mb at high energies, sj 2(2 
Bev)*]. The observed diffraction peak in high- 
energy 7-N scattering suggests that we put in (1) 
dey tap'81?) 4% tap'*1’ Cae 
dQ aQ 





(3) 


g(t) must fall off faster than t™ at large space- 
like ¢ so that 4 jap(s1) remains constant as ob- 
served’ in the several-Bev region, =5-8 mb. We 


choose 
att) -1/| ( 70 me (*) 


for simplicity. The results are insensitive to 
the specific form assumed here.” Inserting (2), 
(3), and (4) into (1) and approximating |s, /t| >1, 
since the second factor (2) weights heavily to high 
scattering energies, we obtain 


x, eto 1S). 





A 


This strong weighting of the amplitude at A to 
large energies s, is what leads to a bump as ob- 
served in the CERN experiments. Of the nine 
final momentum components P , a and k, only 
five are kinematically independent. If we observe 
a final nucleon at fixed Pf this leaves a two-dimen- 
sional integral, 


3,73b a4 -h -a)= ld2ad*h 84 oh-« 
[d'qah Olq+h+p, p; q,) Jd°ad®k6 (k+d, p; 4), 


to be carried out. What (5) tells us is that the 
contribution from large values of s, = (pp +k)? 
=(p;+A)? is greatly enhanced in the phase-space 
integrals. This means that Pi and A tend to be 
antiparallel, and the dominant contribution to 
the observed inelastic scattering comes when 
the target nucleon is hit head-on in the lab sys- 
tem. If the angles of A are constrained in this 
way, we are then led to a unique relation between 
the magnitude of the momentum transfer A?« -/A|’ 
for |A4?|<M?, and the final nucleon energy Ey, oF 


inelasticity,® € 
1Al t( 1 1) 
ie 6 
2 2\IAl m (6) 


Integrating over the range of possible [Al for a 
given E; and 6, leads to a corresponding range of 
preferred € values according to (6). For the con- 
ditions of the CERN experiments these preferred 
inelasticities occur at } Bev<¢<2 Bev and in the 
range of the observed quasi-elastic peaks. 

The more accurate results of detailed calcula- 
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tions reproduce these general features. Inserting (5) into the calculation of the differential cross 


section for Fig. 1 gives 


d*g 72 FF E 1 
 — o, .°(s we 
dE dQ 8(2n)*u? tot 1'°°’E Ip. -p| 
f f t “it “f 
where 
2m? -(m +€)t ~Iz .% | 2me 
— 


and gis the azimuthal angle of @ relative to the 
plane of pj and G;, which determines the polar 
axis. The integral dy is best carried out in the 
system in which the pion and recoil nucleon come 
off back to back; i.e., k +q= 0. The form factor 

F(A?) is introduced at each vertex to cut off the 
integral d(-A*) at maximum values of -A? = (3u 
to 4u)? compatible with the peripheral collision 
assumption of one-pion exchange with a physical 
amplitude inserted at vertex A. The observed 
quasi-elastic bump comes about because the 
emerging nucleon wants to steal most of the en- 
ergy as indicated by (6). On the other hand, the 
integration interval in (7) vanishes when b <e ~—0 
and this leads to a very sharp dropoff of the cross 
section at the inelastic threshold ¢€ =0. 

The second characteristic feature, of a decrease 
in peak height with increasing momentum transfer, 
is assured both by the diffraction factor g(t) in (7) 
and by the integration limits, (8). As -t increases, 
so does a and, therefore, (-A*)min; this reduces 
the integral d(-A”) and further decreases d’o/ 
dEfdQy. The first characteristic feature, that 
the quasi-elastic peak is separated from the elas- 
tic one by a constant € independent of momentum 
transfer ¢ in the range 3 (Bev)?< -t<2 (Bev)?, is 
also established. In order to understand this re- 
sult, note that, according to (8), (-4*) nin in 
creases with -i; this leads to a corresponding 
increase in the mean value of |A| contributing 
to (7). As suggested by (6), an appropriate in- 
crease of |A| with -¢ can lead to a constant € 
independent of ¢. This is what we find to occur 
to a good approximation. 

Calculated curves are shown in Fig. 3 together 
with the data. The calculated numbers correspond 
to a cutoff in the integral Jd(-A?) at a maximum’°® 
of (4u)*. Were we to run the integrals out to the 
limits with F(A?) 1, the characteristic bump 
would remain but its height increases by a fac- 
tor =6 for the 40-mr points and for the calcula- 
tions at E = 16.0 Bev and 6,=56 mr. For the 





a+b (-a2) 


2a 
; d( “09 A aetian dy(s, -m’)’, (7) 





points with the largest momentum transfers, at 
E;=26.10 Bev and 6,=56 mr where -t=2 (Bev), 
the calculations are much less reliable. We are 
farther from the one-pion exchange pole, since, 
for €=1 Bev, |Amin! =4u and the calculated 
cross sections are cutoff-sensitive and reduced 
by an order of magnitude. The neglected proces- 
ses of Fig. 2 may then be relatively important 
and not negligible. 

From the above arguments it appears that the 
process of Fig. 1 is the dominant peripheral con- 
tribution to quasi-elastic small-angle nucleon- 
nucleon scattering at high energies. The contri- 
bution of Fig. 2(b), which was discussed as a 
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FIG. 3. Comparison of observed! cross sections in 
the region of the quasi-elastic bump for different values 
of the experimental parameters with calculations from 
Eq. (7). 
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test of the peripheral approximations,‘ is less 
important for the very small inelasticities, €, 
considered here. This conclusion is also true 
as 6,~0 and we approach the one-pion exchange 
pole.” It wins out only for larger € which keep 
the amplitude of Fig. 1 from developing its dif- 
fraction peak. In an interesting letter which 
helped motivate this work, Selleri* has shown 
that at much lower energies both processes are 
important. 

Our model offers no explanation of the camel’s 
hump structure suggested in the recent experi- 
ments in progress at CERN.’ Its existence, if 
confirmed, may indicate final-state interactions 
or other effects not included in the simple pic- 
ture. 

Finally we note that similar considerations ap- 
ply to high-energy 7-N scattering cross sections. 
For the same kinematics as in the CERN meas- 
urements, there should be a quasi-elastic dif- 
fraction peak in the process 1+N —1+7+N, whose 
height measures the square of the total 7-7 cross 
section in place of the 7-N cross section o¢,¢7(s}); 
in (7). 
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the Air Force Office of Scientific Research. 
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K, DECAY AS A TEST OF UNIVERSAL V-A LEPTON COUPLING* 


M. Bolsterlil and D. A. Geffen 
School of Physics, University of Minnesota, Minneapolis, Minnesota 
(Received July 28, 1961) 


3 


In this note, we show that the experimental 
spectrum of electrons from Kg3 decay provides 
a direct test of universal vector -axial vector lep- 
ton coupling. The test is independent of the K-17 
form factor. We also discuss briefly the best way 





of treating the data so as to obtain information 
about the form factor itself. 

Ke3 decay has been treated extensively in the 
literature.' With the usual V-A lepton current, 
the matrix element for the decay is 


T = (20)*(G/V2)6*(K - p-k = q)(4K oo) “UR Ha (y+ K + a), (Q?) + 30K - 7 af_(Q”)}(1 + y5)0(d). (1) 


We use units in which h=c=m,=1. The mass 
and four-momentum of K, 7, e, v are (M,K), 
(1,q), (meg,p), and (0,k), respectively; 


Q?= (K-q)?= (K-q),?- (K-4); (2) 
f, and f_, the K-7 form factors, are analytic ex- 
cept for a branch on the real axis for Q? > (M+1)?; 


they are real in the range of Q? which occurs in 
Ke3 decay: (M?-1)/2M<Q<M. The equations 


ie)y-k=0, y+ pu(p)=m,v(p) (3) 


can be used to write Eq. (1) in the form 


T = (21)*(G/V2) (4K go) *f,.(Q7)5*(K -p-q-k) 
x u(R)y-K (1+ y,)u(p). (4) 


In Eq. (4) and in the following, we set mg=0. This 
neglect of me is justified as long as |f_/f,| 
«M/me ~10°; the near equality of the Kg3 and 
K ,3 rates’ shows that this latter condition is 
satisfied. 

We find for the number of electrons per unit 
time in the energy interval dp emitted by K’s de- 
Caying at rest 


N(p)dp = dp[G?/2(2a)8|M(M-2p)F (x), _ (5) 


where 





2p(p. -p) M? ~1 
omni. a 
a M-2p ” Prax 2M” (6) 
F (x)= {dy -yf?(2My). (7) 


This energy distribution for the electrons occurs 
(with different normalization) for any combination 
of vector and axial vector in the lepton current. 

Note that p is a double-valued function of x. 
Therefore, N(p)/(M-2p)< F(x) must have the 









same value for the two values of p correspond- 
ing to each value of x. This is, therefore, a sim- 
ple condition that the electron spectrum must sat- 
isfy if leptons are universally coupled via a vec- 
tor-axial vector current. This condition is inde- 
pendent of the form factor which occurs in F(x) 
and depends only on the assumption of V-A lep- 
ton coupling. The data available at present are 
not sufficiently accurate to allow any conclusion 
as to whether or not this condition is satisfied. 
F(x) also has the properties 


F(x)> 0, 
F'(x)>0, F’(e)2=x"F(x), 
F"(x)=f2(2Mx). (8) 


From these properties, it follows that N(p) can- 
not have minima or maxima unless 


p=0, p=p i OF (M?-1)/4M <p < (M-1)/2. (9) 


Thus, except for the end points, maxima and min- 
ima must lie in the range 


114 Mev< p< 180 Mev. (10) 


The histogram given by Furuichi et al.* does not 
appear to satisfy this last condition, although im- 
proved statistics are needed before a definite con- 
clusion can be drawn. 

With sufficient data, it will be possible to deter- 
mine F(x) <N(p)/(M-2p). Besides the single- 
valuedness of F(x), the conditions (7) can be veri- 
fied and, in principle, £?(2Mx) can be determined 
from the curvature of F(x). Note that F(x) can be 
determined more accurately than N(p), since each 
value of F(x) comes from both values of p corres- 
ponding to that x and therefore has an accuracy 
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greater than either of the corresponding values 
of N(p). 
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POSSIBLE EXPLANATION OF THE HIGHER PION-NUCLEON AND K~-p RESONANCES 
IN TERMS OF INELASTIC THRESHOLDS 


James S. Ball and William R. Frazer 
University of California, San Diego, La Jolla, California 
(Received July 31, 1961) 


We have found by means of partial-wave disper - 
sion relations that a large, rapidly rising inelas- 
tic cross section, such as is observed in the vi- 
cinity of the second and third pion-nucleon reso- 
nances,’ can account for a sharp peak in the elas- 
tic scattering. We have calculated the inelastic 
contribution to the higher partial waves by means 
of the strip approximation to the Mandelstam rep- 
resentation,”** in which the principal mechanism 
is the production of the J=1, 7=1 pion-pion reso- 
nance. Assigning to this resonance a position and 
width in accord with recent experimental evi- 
dence,* we find strong inelastic scattering in an- 
gular momentum and isotopic spin states which 
phenomenological analyses have suggested as ac- 
counting for the resonances.*® Moreover, reso- 
nances seem to be implied in approximately the 
correct positions, although we are unable at this 
stage to make really quantitative calculations. 
The same strip approximation calculation seems 
to imply a K~-p resonance, primarily in the D,, 
state, at the observed energy,’ provided the re- 
cently observed K* has J=1 and J=3.° 

In order to find the effect of strong, rapidly 
rising inelastic scattering in a given partial wave, 
consider the amplitude f(v), which assumes the 
following form in the physical region: 


flv) = (e721) 23. (1) 


Above the inelastic threshold v;, 6=5p+i5,, with 
5;> 0 according to unitarity. The variable v is 
chosen appropriate to each particular reaction; 
for equal-mass scattering v =k’, the square of 
the momentum in the center-of-mass system. 
Note that 


Im f = |f|? + (1-7) /4, (2) 
where n=e ~25/, Thus the inelastic contribution 
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to the cross section is proportional to (1 -7”)/4, 
so that 67 is determined if the inelastic scatter- 
ing in the partial wave is known. Let us assume 
for the moment that this is the case. 

In order to keep the argument from being ob- 
scured by inessential complications, we shall 
now make two simplifying assumptions which 
could be relaxed without difficulty. We assume 
that A(v) =f(v)/k is the function which is analytic 
in the v plane except for the singularities given 
by the Mandelstam representation, thereby ne- 
glecting some relativistic effects. Since the point 
we wish to make concerns the effect of the inelas- 
tic scattering on the elastic at rather high ener- 
gies, we also neglect the unphysical cuts. 

We now have the following mathematical prob- 
lem: Find a function A =f/k with the properties: 
(a) It is analytic in the entire v plane, except for 
a cut along the positive real axis. (b) It is real 
on the negative real axis. (c) The function f(v) 
assumes the form of Eq. (1), with given 67, when 
the positive real axis is approached from above. 

The solution is remarkably simple; construct 


oo 5 (v’) 
o(v) = f pw 
Vv. 


, k'(v' =v) (3) 
Using this formula to define 65 in the entire v 
plane, we assert that Eq. (1) provides a solution 
to the problem posed.*"° It can easily be verified 
that A(v) has only the desired branch cuts, pro- 
vided k is defined to have its cut along the posi- 
tive real axis. 

Let us now consider the effect upon the elastic 
amplitude of a rapidly rising inelastic cross sec- 
tion which rises to a value near its unitarity lim- 
it, such as is shown by the solid line in Fig. 1. 
Equation (3) implies that 5p will be a sharply 
peaked function in the neighborhood of the rise 
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FIG, 1. A fit to the Dy, pion-nucleon resonance. The 
dashed line is the calculated o;(Dy,). The dash-dot line 
is the absolute square of the elastic amplitude calculated 
from the inelastic cross section shown by the solid line. 
This inelastic cross section, which joins smoothly with 
our calculated values, has been chosen arbitrarily to 
produce a fit to the observed elastic scattering.’ The 
dotted line is the unitarity limit for the absolute square 
of the elastic amplitude. The variable w is defined as 
w=W-m, where W is the c.m. total energy and m is 
the nucleon mass. 


in the inelastic cross section. The resulting 
elastic amplitude for a particular choice of 6, 
is shown in Fig. 1. Note that the trailing edge 
of the peak is controlled mainly by the unitarity 
limit. The height and exact shape of the peak de- 
pend not only upon the rate of rise of the inelastic 
cross section, but also on its behavior at higher 
energies. It is possible, as illustrated in Fig. 1, 
to find reasonable inelastic cross sections which 
imply elastic peaks sharp enough to fit the higher 
resonances in pion-nucleon scattering." The in- 
elastic cross section itself need not be sharply 
peaked to produce a sharp, sizeable peak in the 
elastic amplitude. This statement is in principle 
subject to experimental verification; for example, 
by measuring the absorption in the Dy, 1=%, pi- 
on-nucleon state above the second resonance. 
Note that our results do not imply a large elastic 
peak at all inelastic thresholds; the condition of 
a rapid rise to near-total absorption must be sat-~ 
isfied. 

We turn now to the second part of the problem; 
namely, to calculate the inelastic cross section 
in each pion-nucleon partial wave and thereby ob- 
tain a value for 6; to substitute in Eq. (3). Asa 
first approximation to the problem, we use the 
“strip approximation” to the Mandelstam repre- 
sentation, in which only those parts of the double 
Spectral functions corresponding to the lowest 
values of the momentum transfer are calculated.*** 














FIG, 2. Diagram giving rise to the pion-nucleon spec- 
tral function in the strip approximation. 


The formulas, which are obtained by applying uni- 
tarity to the Mandelstam representation,’ corre- 
spond to the process shown in Fig. 2. For the 7-7 
amplitude we use a J=1, J=1 resonance formu- 
la.°4* Having calculated the double spectral 
functions for low momentum transfer, we sub- 
stitute them into the equation for the imaginary 
part of the scattering amplitude, 


ImA(s, ¢) =2 [ art's: t) (4) 





t'-t’ 


where the lower limit ¢, is given by the Mandel- 
stam boundary curve. We calculate the spectral 
function a(s, t) correctly out to the boundary curve 
corresponding to four-pion exchange. The effect 
of neglecting these higher-mass-exchange contri- 
butions will be small for a sufficiently high partial 
wave. The J=3 states will certainly be unreliable, 
since the integral in Eq. (4) diverges logarithmic - 
ally if we approximate the pion-pion amplitude by 
the P wave alone. For these states the contribu- 
tions we have neglected, such as the contribution 
of the (3,3) resonance, will surely be important. 
Therefore our calculation has some chance of re- 
liability only for J>3. We project out partial 
waves from Eq. (4) and equate our calculated val- 
ue of Im/; with (1 -7?)/4, since the scattering via 
intermediate two-pion state is inelastic. Let us 
introduce the notation o7(D,,), for example, for 
the quantity (1 -7)/4 in the D,, state. Note that 
this quantity differs from the actual inelastic 
cross section by a factor 47(J +4) /k?. 

Our results for these states are as follows: 
o7(D,,.) rises rapidly in the region of the “thresh- 
old” for production of a 7-7 resonance, but rises 
to a height exceeding the unitarity limit by almost 
an order of magnitude. Our approximation has 
broken down because we have not taken unitarity 
into account for the processes 1+N -7+17+N and 
m+m1+N—17+17+N. It may be possible to remedy 
this defect by means of an approximate treatment 
of these processes in which the particles are al- 
lowed to interact only in pairs. Such a calculation 
is now in progress. If we make the plausible con- 
jecture that in a correctly unitarized calculation 
our large o7(D,,) will remain sizeable and will con- 
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tinue to rise rapidly, we can make a rough calcu- 
lation of the position of a peak in the elastic D,, 
amplitude. With a 7-7 resonance energy around 
5m, this peak coincides with the second reso- 
nance. Moreover, the effect of unitarity should 
not affect the results of our present calculation 
greatly as long as the result is well below the 
unitarity limit. Hence we can calculate the in- 
elastic cross section in the region where it is 
beginning to rise and therefore, according to 
the arguments above, we can roughly predict 
the steepness of the trailing edge of the elastic 
peak. In Fig. 3 these results are compared with 
experiment. Also shown is the inelastic contri- 
bution to the only other state with J > 3 which is 
sizeable in our model, the F,, state. It rises 
rapidly at an energy coinciding with the third 
resonance. These angular momentum assign- 
ments seem to be consistent with experiment.® 
The way they arise in our model can be under- 
stood quite easily: Just above the “threshold” 
for production of a 7-7 resonance, the state in 
which the nucleon and the 7-7 resonance are in 
a relative S wave dominates. This state couples 
with the pion-nucleon system in the Dy, and S, 
states. We are unable to calculate o;(S,,.) relia- 
bly, but we do expect o/(D,.) to rise rapidly in 
this region.’® At a higher energy (thus explain- 
ing the separation between the two resonances), 
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Onn (T= elastic) in mb 
+ 
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20 + 














3.0 70 
FIG, 3. The second and third pion-nucleon resonances. 

The dashed lines are the unitarity limits for elastic scat- 

tering in the Dy. and Fg» states calculated from the dou- 

ble spectral functions. The dotted line is the unitarity 

limit if a background of a saturated 07(Dy,) is added. 

The points are from Falk-Vairant and Valladas.' 
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the relative P wave of the nucleon and 7-7 sys- 
tem will become important. This couples to the 
pion-nucleon P,,., Py, and F,, states. Again we 
cannot calculate o7(P,.). We find that o7(F,.) does 
indeed become important in our calculation, but 
of(Py2) does not seem to attain a sufficiently high 
value to produce a large effect on the elastic scat- 
tering. This numerical result could be changed 
by a more definitive calculation. 

The isotopic spin dependence of the diagram in 
Fig. 2 is such that the J= 4 state is favored over 
the J = 3 by a factor of four to one. This factor 
reduces ok y2) sufficiently that no analog of the 
third resonance is expected in 7*-p scattering. 
However, o7Dgp) is still large enough that the 
elastic amplitude in this state should show a peak 
around 800 Mev in 7*-p scattering. This may be 
the knee in the cross section which Carruthers 
has identified with the D,,, state.® 

We have made a similar calculation for K~-p 
scattering around 1 Bev/c lab momentum, where 
a peak in the cross section has recently been 
found.” Again we calculate the diagram in Fig. 2, 
but now the incident and final mesons are K~, and 
the intermediate resonance is the K* resonance 
found by Alston et al.**!® If we assume that the K* 
has J=1 (experiments have not yet distinguished 
between J=0 and J=1), then we again find a 
strong, rapidly rising o7(D,,.) at just the right 
position to correspond to the observed peak (see 
Fig. 4). This implies a nearly saturated D,, con- 
tribution to K* production at 1.15 Bev/c, the ener- 
gy in the experiment of Alston et al, Such a con- 
tribution to the process K~ +p~—K°+7~ +p turns 
out to be only 1.0 mb, however, which is not in 
disagreement with the experimental value which 
is of the order of 1 mb. In this process o7(Fy,) 
does not come out to be strong, largely because 
of the narrowness of the K* resonance as com- 
pared to the 7-7 resonance. We see from Fig. 4 
that with a reasonable background the D,, state 
does not seem to be able to account for the entire 
peak. The difference is probably due to the S,, 
state which we are unable to calculate reliably. 

If we assume the K* has J=0, we find no large 
amplitudes with J>3, because a nucleon and K* 
in a relative S state can then couple only to J=3 
(which is clearly ruled out by the height of the 
observed peak). The relative P state does not 
contribute strongly for this process. Thus we 
definitely require the K* to have J=1 if our model 
is to be valid. If we further assume /= 3, as in- 
dicated by the K* experiment, we find that in K-N 
scattering the J=0 state is favored over J=1 bya 
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FIG. 4. The K~-p total cross section.’ The solid 
line is the calculated unitarity limit for elastic scatter- 
ing in the Dy. state, plus a constant background of 30 
mb. The dashed line is a “fit” to the data, obtained by 
a suitable continuation of the calculated inelastic cross 
section. 


ratio of 9/1. This agrees with the fact that no 
large peak is seen in the K~-m cross section. Our 
model implies that there should also be a peak in 
the K*-n cross section, but it is impossible to 
predict its height from the present crude calcula- 
tion. It is probably considerably smaller than the 
K~-p peak, because there are more inelastic chan- 
nels open in the latter process, which may help to 
raise o7(D,,) close to the unitarity limit. 
Additional resonances are probably implied in 
1™-Y and K-Y systems, although detailed calcula- 
tions must be carried out, which will depend upon 
assumptions about strange-particle couplings. 
The positions of these resonances will be around 
center-of-mass total energy mz-,+my and mx* 
+my, respectively. 
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LEPTON PAIRINGS IN THE TWO-NEUTRINO THEORY 


G. Feinberg,* F. Gtirsey, f and A. Paist 
CERN, Geneva, Switzerland 
(Received July 31, 1961) 


In this note we will discuss some possible impli- 
cations of the assumption that two distinct neu- 
trinos are involved in the weak interactions. In 
particular, we consider the possibility that the 
association of the neutrinos with charged leptons 
in strangeness-changing decays is different from 
the pairing in strangeness-conserving decays. We 
will see that this possibility, if fulfilled, would 
have important consequences for proposed inverse 
neutrino reactions. 

We recall that the existence of two neutrinos 
was assumed in order to forbid the decay u+e+y 
by a selection rule,’ according to which the » and 
e are assigned different values of some quantum 
number. In order to make this consistent with the 
ordinary » decay, two distinct neutrinos v, and v, 
with the same helicity are introduced, and the yu 
decay is assumed to be of the form 


pe" +, + Vy (1) 


This follows in a natural way from the common 
assumption that all strangeness-conserving decays 
arise from an interaction of the form J1J » pro- 
vided that J=np...+e"v,+ u~v_. Such a structure 
implies also that the 7 decays are 


1 +et+v,, (2) 
1* » t+ Vo. (3) 


These processes will not lead to 1 +e+¥y to any 
order.’ 

It has been usually assumed that the pairings 
(uv,) and (ev,) are, in fact, universal and also 
occur in strangeness-changing reactions. The 
purpose of this note is to point out that if two neu- 
trinos exist at all, there is the alternative pos- 
sibility that in decays with AS=+1, the pairings 
(ev,) and (uv,) may occur. In such a theory the 
two-body leptonic decays of the K* would be 


Kt + ur+ Vi» (4) 


set+ vy. (5) 
We will refer to such decays as “neutrino flip” 
decays. In the conventional theory the decays 
would instead be 


Kt > pt+ V2» (6) 
»et+ V}- (7) 


Clearly, no direct experimental evidence exists 
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at present to decide between the reactions (4), (5) 
and (6), (7), although the arguments on forbidding 


u.+e+y to all orders implies that they cannot co- | 


exist. If the reactions (4), (5) rather than (6), (7) 
occur, there will be important experimental con- 
sequences, some of which we indicate here. 

(1) If high-energy neutrinos are obtained from 
the decay of a beam of positive pions, which has 
a “contamination” of positive K mesons, the net 
neutrino beam would be a mixture of v, and p, 
and therefore could generate both the reactions 


v»tnre-+p, (8) 


and 
Vatn> +p. (9) 


We have been informed? that in neutrino experi- 
ments in progress, the K contamination of the pion 
beam may be ~20%. In the presence of neutrino 
flip decays this might lead to a sizable electron 
production. If an experiment is performed with 
neutrinos from a contaminated beam and no elec- 
trons are found, it would not only establish the 
two-neutrino theory, but also rule out the exist- 
ence of neutrino flip K decays. Conversely, the 
mere observation of electrons produced by neu- 
trinos would disprove the two-neutrino theory 
only if a “pure” pion beam was the source of neu- 
trinos. (We neglect the one neutrino in 10* coming 
from 7 +e+v decay.) 

(2) A related point is the availability of high- 
energy v,’s from K* decay if the neutrino flip de- 
cay (4) occurs. Such high-energy v,’s are not ob- 
viously available from any other source and would 
be a good tool in studying the symmetry between 
muon and electron interactions, by comparing for 
example the reactions (8) and (9). It will there- 
fore be useful in this connection to be able to dis- 
tinguish experimentally between inverse neutrino 
events induced by K* neutrinos and by 1* neutrinos. 
For this purpose kinematic discriminations in a 
“contaminated” beam may be sufficient. 

(3) The occurrence of neutrino flip K decays 
necessarily implies the occurrence of neutrino 
flip in hyperon-lepton interactions with AS= AQ. 
This means that in a neutrino absorption experi- 
ment, where the neutrinos are known to be of one 
type, the charged leptons produced with hyperons 
will not be the same as those produced with nu- 
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cleons. In particular the reactions 


Va+p>A°+e*, (10) 

V,+p>A°+ pt, (11) 
will occur, rather than 

Vat p~A°+ ut, ete. (12) 


In addition, we would like to make some com- 
ments of a more theoretical nature. 

(a) If the lepton pairings in strangeness-chang- 
ing decays are indeed different than those in 
strangeness-conserving decays, it would suggest 
that leptons carry a quantum number somehow 
related to strangeness,* and that the weak leptonic 
interactions are characterized by a selection rule 
involving both of these numbers. In this connec- 
tion we mention that the reactions (4), (5) also im- 
ply the existence® of p*+e*+v,+7,. These points 
will be taken up elsewhere. 

(b) Unless special parity structures of the lep- 
tonic, strangeness violating, couplings prevail, 
it follows that if the Kj2 modes (4) and (5) do show 
neutrino flip, all AS/AQ=+1 transitions do the 
same. Thus in the K;3* decays and in the (|AS|=1) 
B and » decays of A, =", and =” we would expect 
also (e,v,) and (u,v,) pairings. The same is true 
inK°+m~+et+y,. If AS/AQ= -1 transitions occur, 
it is evidently impossible to determine the nature 
of the lepton pairing in such reactions from Eqs. 
(4) and (5). To reconcile the absence® of nonlep- 
tonic AS= 2 transitions with a relatively simple 
current-current structure, it is attractive to con- 
sider the possibility that, if reactions with AS/AQ 
=~1 take place at all, the pairing is (e,v,), (u, v,). 
In this case we would have 


K° +1" +et+v, (13) 
>m*+e7+0,, (14) 
R°+1*+e7+0, (15) 
+1" +eT+y,. (16) 


We now assume CP invariance. Then regardless 
of the relative pairings in the reactions (14), (16) 
as compared to (13), (15), we have 


(17) 
(18) 


R,(a)=R,(a*), 
R,(0 *) ” R,(n*), 


where R, (R,) denote the total rate of decay of K, 
(K,) into the 7 with charge as indicated plus the 
corresponding electron-neutrino pair (and like- 
wise for the ,.-neutrino pair). We now consider 
the following three instances. 


(i) No AS/AQ= -1 transitions. We have now in 
addition 


R,(7~)=R, (a7). (19) 


(ii) AS/AQ= -1 transitions occur and the same 
lepton pairings (e,v,), (u,v,) appear in both 
AS/4Q=+1. This means that in Eqs. (14) and (16) 
v, and v, occur, respectively, rather than y, and 
v,. The relation (19) is now in general not true 
due to different coherence effects in K, as com- 
pared to K, decays. 

(iii) AS/AQ = -1 transitions occur and the lepton 
pairing is as indicated in Eqs. (13)-(16). Now, for 
example, the K, decay into 7~ +e* + neutrino is the 
incoherent superposition of the channels (13) and 
(16). Thus the rate is one half the sum of the par- 
tial rates. As a result the relation (19) is true 
even in the presence of (CP-invariant) AS/AQ =+1 
transitions.’ Thus we have here the possibility of 
an experimental distinction between the cases (ii) 
and (iii). Observe also that if the AS/AQ =+1 tran- 
sitions satisfy more specifically |A7T|=4, then the 
following conclusion is immediate: If AS/AQ= -1 
transitions occur, with pairings as in Eqs. (14) 
and (16), then the ratio of rates (K,°+1*+e* + v)/ 
(K* + 7°+e*+ v) is larger than two. 

It should be noted that these last results depend 
only on the lepton pairing in AS/AQ =+1 relative 
to AS/AQ=-1. It is independent of the lepton pair- 
ing in AS/AQ=+1 relative to strangeness-conserv- 
ing reactions. 

In addition it should be emphasized that if AS/AQ 
=-1 involves (e, v,) and (u,v), as in Eqs. (14) and 
(16), then for example v,+n+2*+e™ remains for- 
bidden, while for pairing (e, v,), (u,v,) this reac- 
tion is allowed. 

(c) It has been surmised that all the weak inter- 
actions are generated by a sum of product of cur- 
rents such as ($+j+s)T-(J+j+s)+(jo+ s,)t- (jot So)- 
This extension of the universal Fermi interaction 
to include strangeness-changing processes has 
sometimes been questioned, because of the ab- 
sence of neutral lepton currents and because of 
the slowness of hyperon decay into leptons com- 
pared to the universal rates.* Regardless of these 
problems we would like to emphasize that the ex- 
istence of the neutrino flip K decays would rule 
out this particular current-current structure, since 
different lepton currents occur in AS=0 decays 
than in |AS|=1 decays. Of course, a more com- 
plicated current structure can be envisaged. 

(4) If all weak interactions including nonleptonic 
ones are considered to all orders of perturbation 
theory, then the absence of u +e+¥y decay rules 
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out neutrino flip K decays.® This follows because 
the K meson can convert into a pion through virtual 
weak nonleptonic interactions, and therefore the 
chain p~K+v,+17+v,+e+y can occur to third 
order in weak interactions. From this standpoint 
there is therefore a cogent objection to our pro- 
posal, if the cutoff of weak interactions is very 
large. However, we believe that this theoretical 
question is sufficiently open to warrant a direct 
experimental test of the existence of neutrino flip. 

We are grateful to CERN for its kind hospitality 
and support. 

After this note was submitted, Dr. B. d’Espagnat 
kindly brought to our attention a preprint’® by Dr. 
S. Bludman in which the possibility of neutrino flip 
for a AS= AQ is already considered. 
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The ng or Am resonance at 1385 Mev has been 
observed by several experimental groups.’ One 
of the theoretical explanations of the Y Fg is that 
it is a quasi-bound state of the s-wave KN system. 
In fact, Dalitz and Tuan,? from their rough, con- 
stant-scattering-length analysis of the s-wave 
Kp data, predicted that such a resonant state 
might exist. Discovery of the Y,* led to a discus- 
sion by Ross and Shaw’ of the possibility of putting 
the association of the Y,* with the KN s state toa 
detailed test by means of an effective-range analy- 
sis, and as a result determine the KYN parities. 

The purpose of this note is to report on an ex- 
tensive survey of effective-range solutions of the 
s-wave Kp data. The results can be summarized 
as follows: Assuming that the Y * is associated 
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with K“p s state, the only acceptable solutions 

found were for the D7 channel in a p state (thus 
even KEN parity), whereas the A7 channel could 
be either in an s or p state. In addition, among 
three different sets of effective ranges, only the 
one corresponding to a moderately short “range 
of forces,” i.e., 0.4 f, led to satisfactory solu- 





tions. 

The scattering amplitudes Tj; for a system of 
n coupled two-particle channels in a given state 
of total angular momentum J, z component Jz, 
isotopic spin /, parity P, etc., can be expressed 
as* 


1+1/2  2l+1,-1 


(M -ik )ok 


The partial cross section from channel i to chan- 


T=k +1/2 (1) 
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nel j is given by 


_ 48 8341.6, os 
"ij ke 2s.+1 TNs 





(2) 


where S; is the spin of the initial (unpolarized) 
channel. The (xm) matrix M is an even function 
of all the channel momenta k; and thus remains 
real and symmetric when all the channels are not 
open (below the ith threshold k; ~ix;, x;> 0). The 
diagonal matrix k/ + 1/2 has entries 


git 1/2 lj+1/2 


)..=5.R. ; 
yj yi 

where J; is the angular momentum of the ith chan- 

nel. According to the effective-range theory de- 

veloped in A, M has the following energy depend- 

ence near any fixed energy E, (which we take as 

the K ~p threshold for the present application): 


M=M(E,)+3 Rk? -k7(E,) J, (3) 


where ® is a real, approximately energy -inde- 
pendent and diagonal (even for 1; #1;) matrix whose 
elements are roughly interpretable in terms of an 
interaction range. We have 
_ . 1-21;,, 2 
M,,=M,Eq)+25,R, Lk; 

where C,=1, C,=-3, and R; is a measure of the 
range of forces in the zth channel. 

The s-wave J=0 KN channel is coupled to the 
=a channel whereas the J=1 system consists of 
three coupled two-body channels (KN, =7, Az). 
In addition to the energy dependence of the iso- 
topic spin amplitudes there are two charge-de- 
pendent modifications which are made: correc- 
tions due to the K"p-K°n mass splitting, and 
Coulomb corrections due to the charged K“p chan- 
nel, 

There are eight pieces of data from s-wave 
K“p reactions which are used in the present analy- 
sis: The capture of K~ mesons from s orbitals in 
(K~,p) atoms yields the production ratios® >*1~/ 
D"n*/E°nr°/An®, At an incident K~ laboratory mo- 
mentum of 175 Mev/c there are measurements 
on the elastic and charge-exchange cross sections 
as well as the production in the charged hyperon 
channels 5*~ and =“a*.® At 400 Mev/c, the p- 
wave amplitudes become dominant, and there is 
a very strong restriction’»® on the s-wave elastic 
Scattering: 


2 
-k Eg), 4) 


0, (400 Mev/c) < 22 mb. 


The position (~ 50 Mev below E,), width (full 
width at half maximum <30 Mev), and (27/Am)7;-4 






ratio (<0.2) measurements on the n° have been 
fairly well corroborated.’ We will assume that 
the Y,* has J=4 and P(KY,*N) is even, i.e., the 
At resonance is associated with the s-wave KN 
system. 

The computational procedure we use is the fol- 
lowing: The 11 s-wave KN data (8 above-thresh- 
old K~p measurements and 3 resonance param- 
eters) are to be fit using the effective-range form- 
alism (1)-(4) properly modified to take into ac- 
count the charge-dependent corrections mentioned 
previously. We use 10 parameters: 9 zero-range 
Mij(Eo) and one effective-range R, i.e., we as- 
sume R;=R for alli. For each possible combina- 
tion of P(KZN) and P(KAN) we considered more 
than 10000 sets of zero-range Mjj(EQ). All these 
sets of M;;(Eg) were chosen so that they fit the 
at-rest data® and yield complex scattering lengths 
a;(E9) for the KN channels of the Dalitz-Tuan (a-) 
type (in order to omit all solutions which seem 
very unlikely to lead to the Y >," Thus only some 
of the parameters were chosen arbitrarily [e.g., 
we let the matrix element 

1+1/2 -1,1+1/2 
& “) =1 . ds »u 
evaluated at E, take on 7 values between +1.2 and 
-1.2], whereas others were either determined or 
more narrowly restricted by these conditions im- 
posed at the beginning of the computer program. 

The relationship between the scattering length 

and the M matrix is, for a given isotopic spin, 


“1/ah)=M,+M ‘2 ms 
ign 
° — y2l+1 
where M,; is the 4j minor of (M-ik ) and x 
denotes the KN channel. [Note that for the one- 
channel system, -1/a(k)=M=k coté.| Thus the 
zero-range approximation, M=M(E,), is not 
equivalent to the constant scattering-length ap- 
proximation, a(k)=a; for constant M, a(k) may 
still have considerable energy dependence. For 
each set of Mij(EQ), we carried the amplitudes to 
other energies according to the three different 
prescriptions: the zero-range approximation 
Mi =Mij(E9), R=0.4f, and R=0.7f. The calcu- 
lations required several hours on a C.D.C. 1604 
computer. 

In the present survey, only 8 solutions satisfied 
the 11 restrictions. All eight solutions were for 
the =7 channel in a p state, i.e., even KIN parity 
[in spite of the fact that twice as many sets of val- 
ues were investigated for odd P(K=N) as even]. 
The number of solutions were evenly divided be- 


(5) 
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tween even and odd P(KAN). In addition, all these 
solutions had a range of interaction R= 0.4 f (yield- 
ing strong energy dependence in the D7 channel). 
Some additional aspects that these “acceptable” 
solutions have in common are as follows: (i) No 
isotopic spin zero resonance was found. The en- 
ergy dependence of 1/a,(k) was very rapid (and 
nonquadratic) and even those solutions which in 
the constant-scattering-length approximation 
would yield a broad J=0 resonance had quite small 
nonresonant J=0 amplitudes below E,. (The pos- 
sibility of having the s-wave KN amplitudes accom- 
modating both the Y,* and a narrow Y,* seems 
completely ruled out.) (ii) The 5~/5* ratio 
changes rapidly from 2.18 at rest to =1 at 175- 
Mev/c incident K~ laboratory momentum, then 
remains fairly constant (slightly more =*n™ than 
Dat). (iii) The A/(A+2°) ratio increases from 
at-rest values =~ 0.2 to ~0.6 at 300 Mev/c. (iv) 
(=/A)z=1, which is <0.2 in the resonance energy 
region and ~1 at threshold, decreases to a value 
«1 at 300 Mev/c.® (v) The charge-exchange scat- 
tering remains fairly large, i.e., it does not fall 
off greatly from its value of ~12 mb at 175 Mev/c. 
With the existing above-threshold s-wave K~p 
data, we find (nonacceptable) effective-range solu- 
tions which fit these data and yield scattering 
lengths near threshold of the Dalitz-Tuan (a-) 
type (which would in the constant-scattering- 
length approximation predict a Am resonance hav- 
ing a position and width corresponding to the ob- 
served Y,*) and yet show no resonant behavior. 
Thus definite evidence for the (a-) type Kp scat- 
tering length solution is not in itself conclusive 
evidence that the Y,* must be identified with the 
s-wave KN system. Additional K~p (and K°p) 
measurements shedding light on the detailed en- 
ergy dependence of the amplitudes, such as the 
0e)(400 Mev/c) datum, which cannot be fit in a 
constant-scattering-length approximation, may 
change this situation. However, identification of 
the Y,* with the s-wave KN system together with 
the present effective-range survey indicates that 
the KIN parity is even. Determination of the an- 
gular momentum of the Az channel, for example, 
from measuring the polarization of the A produced 
in the eg decay would then give the KAN parity 
(there is some indication that the A7 is in ans 
state’). 








The details of this calculation (including a dis- 
cussion of the accuracy and applicability of the 
effective-range formalism for the present prob- 
lem,’° e.g., the restriction |k;?-k;7(Eo)|R?<1 is 
satisfied over the entire energy range considered) 
along with a treatment of global symmetry as ap- 
plied to KN reactions will be published elsewhere, 

The encouraging results of the present survey 
(e.g., only a small number of acceptable solutions 
were found) indicate that a gradient search" for 
the M;j(Eg) and Rj would be warranted in the near 
future. 

We would like to thnak Professor R. H. Dalitz 
for a very helpful discussion. 
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ANOMALOUS THRESHOLDS AND FINAL-STATE INTERACTIONS IN PRODUCTION PROCESSES* 


R. F. Sawyer 
University of Wisconsin, Madison, Wisconsin 
(Received July 18, 1961) 


Evidence for resonant two-body final-state in- 
teractions in production processes with three- 
particle final states has recently been obtained 
in the form of energy spectrum peaks in numer- 
ous experiments involving strange particles.’ 

It is less obvious how the effects of nonresonant 
scattering in strongly interacting three-body 
systems are to be estimated, and under what 
circumstances they can be expected to provide 
the dominant dependence of a production ampli- 
tude on some variable. Nevertheless nonresonant 
S-wave pion-pion effects are now suspected to be 
the cause of the strange He* momentum distribu- 
tion in the p+d experiment of Abashian, Booth, 
and Crowe.”** If this explanation is correct, a 
similar effect should be observed in the final nu- 
cleon spectrum in the process N+1+N+7 +7. 

It is the purpose of this note to point out some 
complicating effects, due to anomalous thresh- 
olds, in the latter process. These effects will 

be common to many production processes (in- 
cluding p +d + He* +27°) and may muddle the in- 
terpretation considerably. The virtue of the par- 
ticular process, N+7+N+m+17, is that it allows 
a quantitative estimate of the difficulties. 

We shall devote our attention to the dependence 
of the production amplitude on a variable, q, 
which is the momentum of one of the final pions 
in the center-of-mass system of the two outgoing 
pions. Only pairs of pions in the T=0, J=0 state 
and only small values of g will be considered. 
Whatever the complete analytic structure of the 
production amplitude, there should be a cut in 
the q* plane running from zero to infinity, orig- 
inating in the 7-7 rescattering. This cut may be 
taken into account by writing for the production 
amplitude, A (all variables except q’* suppressed),' 


° 1f{~ (9) , 
A(q’) =M(q’) oot f qu -g-tae* . 
The second factor here is the inverse of the stand- 
ard denominator function, D,* computed from the 
T=0, S-wave 1-7 phase shift. M(q*) is analytic 

in g along the physical 7-7 cut. One is now 
tempted to assume (and it is this assumption 

which we shall criticize) that for small g the 
function M(q*) can be approximated by a constant. 
In a scattering-length approximation with an at- 


tractive 7-7 force, the dependence of the cross 
section for small q* would then be as (1+ qa’), 
which for a large scattering length, a, would 
produce an enhancement for small g*. This is 
the basis upon which an attractive, T=0, 7-7 
force can be considered an explanation for the 
T =0 anomaly in p+d collisions.”* 

As has been pointed out,** production ampli- 
tudes may have a wide variety of complex branch 
points. Of these there is one particular kind of 
anomaly which seems to be of importance with 
respect to the above-mentioned arguments on 
final-state interactions. This is associated with 
the reduced graph® of Fig. 1(a). We define varia- 
bles for Fig. 1 as 


We -(k, +p,)’, 
A= -(p, -p,)’, 
s =-(k, +k,)?=4(¢? +”). (1) 


A. For sufficiently small values of W, the dia- 
gram of Fig. 1(a) has only the normal cut in s. 
As W is increased an anomalous branch point 
breaks off from the normal threshold at s =4y? 
and moves to the left. When W becomes unstable 
[W =(M +)*] this branch point moves into the 
complex plane to the point s,,° 


«= E(W+M* ~ yt) 
b 2M" 





Z 
al" -2(M?+yu)W 
+ (M? ~ p?)?}¥9(4M? - ?)¥?, (2) 
B. Let the 7-7 interaction of Fig. 1(a) be local 
(ag*). Then the cut originating at the complex 
branch point, sp, is quite deformable. However 


it must terminate on the normal threshold at s 
=4y? if ordinary unitarity is to be used along the 


Kx 


K, CK --9Ks ia 


rl ye Oe * B 
ve! ‘e OR. 
P P : 


“i: ' 


FIG. 1. Two reduced diagrams for 7+N—7+7+WN, 
Solid lines are nucleons. Dashed lines are pions. 
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normal cut, 4u?<s<~ (this is true at least in 
perturbation theory). The anomalous cut there- 
fore gives a branch point at s =4 which is not re- 
moved by factoring out D™. It remains to esti- 
mate the importance of this singularity. For the 
graph of Fig. 1 the discontinuity across the com- 
plex branch cut may be calculated using gener- 
alized unitarity,® in terms of A, G (1-N coupling 
constant), and the 7-N scattering amplitude. The 
™-N scattering will be on the mass shell, but at 
a complex momentum transfer. We shall return 
to this calculation. 





C. Figure 1(a) seems to be the only contribu- 
tor to the additional singularity at s=4u? men- 
tioned above, though a further check of some 
higher order graphs is in progress. Of reduced 
diagrams which contain the 27 state in the s chan- 
nel, the next to be considered is that of Fig. 1(b). 
This graph is not anomalous in s for any value of 
W. This is by virtue of the internal instability at 
the vertex, B.°® 

D. If all lines in Fig. 1(a) were boson lines, 
the discontinuity across the anomalous cut run- 
ning from s = sy to s=4u? could be calculated 
from the expression,® 


disc A =(KF)AGJd*p 5(p? +M*)d[(p -p,)” + 1? ]6[ (p -Py -k,)? + u*|T(W, - (p, -p)”), (3) 


where T(W, t) is the 7-N scattering amplitude as 
a function of energy and momentum transfer, and 
(KF) stands for kinematical factors. One obtains 


—- T(W, t) 
discA = (KP bony (4) 


where 
p? 2\y/2 2\1/2 

t =-s5[1 -xy ¥ (1 -x?)¥7(1 -y?)*], 

x =(1/aB)(y5 -2sW), 


y =x +(2W/aB)(A - u?), 
a =(y? -4sw)”, 


B = (6? -4y?w), 
y=W-M?* +s, 
6=W -M?+ 7, 


The integral in ¢ is around the cut in the denomi- 
nator. It can be seen that in the high-energy lim- 
it, W>M, A, s, the region in the cosine of the 
m-N scattering angle implied by the integration 
limits in (4) becomes limited to a small region 
near cos@=1, and that the momentum transfer 
becomes real (even for the complex values of s 
under consideration). In fact, for high energies 
this region in angle remains somewhat within the 
forward diffraction peak. Therefore we insert 

in (4) for T(W, ¢) a diffraction approximation for 
a black sphere of radius R, 


iQRJ,[R(-t)*] 


(1 (5) 


T(W, t)= 





E. Considerable complication develops ina 
proper calculation in which nucleon spin and odd- 
parity pions enter. For Fig. 1 alone four ampli- 
tudes must be introduced and the discontinuity of 
each calculated along the anomalous cut. This 
has been done in the high-W limit and in the dif- 
fraction approximation described above. Here 
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we quote only the result. The spin dependence 

is just “(p,)y,u(p,) as in the single-pion pole term. 
The result will be displayed as the ratio of the 
contribution of the anomalous cut to that of the 
pole cerm, thus eliminating A, G, and kinematical 
factors. The isotopic spin dependence of the two 
contributions is identical in our diffraction ap- 
proximation. We obtain 


ve 4 = 
= i(A-p*)(Ru)* s-4p 
= z In . (6) 

8u S-S, 





pole 


For Ru =1 the logarithmic singularity is both 
large and strongly dependent on s near s =4y’. 
This latter property apparently will cause some 
difficulty in using scattering-length approxima- 
tions to treat final-state interactions at low rela- 
tive momentum. 

The domain in which there is such a difficulty 
is better estimated from a somewhat different 
comparison. In the one-pion pole contribution 
to the production process let us insert, instead 
of Ay*, a scattering-length approximation to the 
pion-pion amplitude, 


fag 7/1 -iga). (7) 


The contribution of diagram (la) is calculated as 
before, taking a coupling constant appropriate to 
the scattering length in (7). Now we compute the 
derivatives with respect to q of Agnom and Apole 
in order to give a comparison of the anomalous 9 
dependence with the g dependence of the pole term 
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implied by (7). One obtains 
(4/dq)A. om . (A - u?)(Ru)*(1 -éga)? 


= (8) 
(d/dq)A 





8u*ga 


For 4=0 and Ru =1/2” (as indicated by high- 
energy ™-N total cross sections), we find that 

we must go to ga > 1/32 to see the direct q de- 
pendence from the 7-7 scattering amplitude dom- 
inate the anomalous q dependence. For example, 
for A=10 there is clearly no domain of ga<1 for 
which the ratio (8) is small. The scattering- 
length approximation is never valid for very low 
q. There may or may not be a domain of validity 
for intermediate g, that is, in a region sufficient- 
ly far away from the logarithmic singularity but 
for small enough ga to justify a scattering-length 
approximation. 
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ERRATUM 





SUPERCONDUCTIVITY IN Nb,Sn AT HIGH CUR- 
RENT DENSITY IN A MAGNETIC FIELD OF 88 
kgauss. J. E. Kunzler, E. Buehler, F. S. L. Hsu, 
and J. H. Wernick [Phys. Rev. Letters 6, 89 


(1961) ]}. 


Reference 5 is incorrect. Instead of J. C. Mc- 
Lennan et al., Trans. Roy. Soc. Can. 24, 25 (1930), 
the correct reference should be R. R. Hake, D. H. 
Leslie, and T. G. Berlincourt, Bull. Am. Phys. 
Soc. 4, 362 (1959); de Haas and Voogd, Commun. 


Phys. Lab. Univ. Leiden No. 208b (1930). 





